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ABSTRACT 

The  Geodetic  Absolute  Sequential  Positioning  (GASP)  program,  as  utilized  by  the 
Defense  Mapping  Agency  (DMA),  processes  static  GPS  measurements  collected  with 
the  TI  4100  GPS  receiver  to  estimate  geodetic  point  (absolute)  positions.  In  this  thesis, 
the  GASP  program  is  modified  to  accept  data  from  different  receiver  types,  the  estimated 
point  positions  are  compared  to  positions  produced  by  the  Transit  Doppler  positioning 
system,  the  between-receiver  estimates  are  compared,  the  difference  between  estimates 
using  the  broadcast  and  the  precise  are  examined,  and  the  effects  of  Selective  Avail- 
ability assessed. 

During  the  Monterey  Bay  Precision  Positioning  Experiment  (MBPPE),  conducted 
in  the  Winter  of  1990-91,  a  large  set  of  static  GPS  positioning  data  was  collected  with 
four  types  of  GPS  receiver;  the  TI  4100,  the  Trimble  4000ST,  the  Ashtech  LD  XII,  and 
the  Magnavox  MX4200.  Additional  static  GPS  measurements  were  obtained  with  the 
TI  4100  receiver  at  a  reference  site  established  to  support  the  experiment.  A  third  data 
set  was  collected  after  activation  of  Selective  Availability.  Measurements  collected  with 
the  TI  4100,  Trimble,  and  Ashtech  receivers  were  subsequently  processed  with  GASP 
using  broadcast  and  precise  ephemerides  to  produce  point  position  estimates.  In  order 
for  GASP  to  accept  the  data  from  the  Ashtech  and  Trimble  receivers,  the  program  had 
to  be  modified. 

The  positioning  results  obtained  are  analyzed  for  accuracy  and  precision.  The  ac- 
curacy of  the  GASP  GPS  estimates  is  determined  by  comparison  to  independent  esti- 
mates obtained  by  the  Transit  Doppler  positioning  system.  Precision  or  repeatability 
(i.e.,  consistency  of  the  estimated  positions)  is  also  examined. 

Analysis  of  the  accuracy  and  repeatability  reveals  little  difference  between  the  posi- 
tions computed  for  the  three  receivers  using  the  precise  ephemeris  and  that  all  three 
provide  good  agreement  to  the  Transit  Doppler  positions.  All  three  receivers  are  capable 
of  providing  geodetic-quality  point  positions.  It  is  also  clearly  demonstrated  that  the 
precise  ephemeris  does  produce  a  more  accurate,  higher  precision  solution  than  the 
broadcast  ephemeris.  The  activation  of  Selective  Availability  has  substantially  degraded 
position  solutions  available  from  the  broadcast  ephemeris. 
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THESIS  DISCLAIMER 


The  reader  is  cautioned  that  computer  programs  developed  in  this  research  may  not 
have  been  exercised  for  all  cases  of  interest.  While  every  effort  has  been  made,  within 
the  time  available,  to  ensure  that  the  programs  are  free  of  computational  and  logic  er- 
rors, they  cannot  be  considered  validated.  Any  application  of  these  programs  without 
additional  verification  is  at  the  risk  of  the  user. 
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I.     INTRODUCTION 

A.     THESIS  DESCRIPTION  AND  OBJECTIVES 

1.  Monterey  Bay  Precision  Positioning  Experiment  Overview 

The  Monterey  Bay  Precision  Positioning  Experiment  (MBPPE)  was  designed  to 
assess  the  performance  of  some  commercially  available  Global  Positioning  System 
(GPS)  receivers  and  the  processing  software  developed  to  support  the  receivers.  In 
particular,  the  static  and  dynamic  positioning  solutions  of  the  Trimble  4000ST,  Ashtech 
LD  XII,  Magnavox  MX4200,  and  Texas  Instruments  TI  4100  receivers  were  examined. 
All  but  the  MX4200  are  geodetic  quality  receivers.  The  Trimble,  Ashtech  and 
Magnavox  models  exhibit  some  of  the  latest  developments  currently  available  in  GPS 
receiver  technology.  The  TI  4100  is  an  older  model  but  is  still  widely  used  throughout 
the  geodetic  community. 

The  overall  objective  of  the  experiment  was  to  acquire  a  large  data  set  of  GPS 
measurements  in  both  static  and  kinematic  modes  with  the  different  receivers,  to  process 
the  data  by  utilizing  processing  techniques  and  software  relevant  to  a  given  application, 
to  interpret  the  results  from  a  position  accuracy  and  error  analysis  perspective,  and  to 
evalute  the  receivers  and  processing  techniques  based  on  these  results. 

2.  Thesis  Overview 

The  purpose  of  this  investigation  is  to  determine  whether  meter  level  point  po- 
sition accuracy  is  attainable  with  the  TI  4100,  Trimble,  and  Ashtech  receivers.  The 
Magnavox  receiver  is  not  considered  in  this  study  because  it  is  a  single  frequency  re- 
ceiver, whereas  the  other  receivers  operate  on  dual  frequencies.  This  is  a  serious  limita- 
tion in  determination  of  point  positions  since  the  dual  frequency  correction  cannot  be 
applied  to  correct  for  ionospheric  refraction. 

This  study  is  concerned  with  one  aspect  of  geodetic  positioning,  the  absolute 
determination  of  a  point  position  in  some  commonly  used  reference  frame.  There  are 
two  general  categories  of  positioning,  static  and  kinematic.  As  may  be  inferred,  static 
positioning  involves  determining  the  positions  of  stationary  objects  and  kinematic  posi- 
tioning the  positions  of  moving  objects.  Static  positioning  may  be  divided  into  relative 
and  absolute  positioning  subclasses.  Relative  positioning  is  described  as  the  estimation 
of  the  vector  (baseline)  connecting  a  known  station  to  an  unknown  station.  Absolute 
or  point  positioning  may  be  described  as  the  estimation  of  the  vector  connecting  the 


origin  of  a  global  reference  frame  with  an  unknown  station  [Ref.  1:  p.  1],  It  then  follows 
that  relative  kinematic  positioning  may  be  described  as  the  estimation  of  a  vector  con- 
necting a  known  static  station  to  an  unknown  moving  object.  Station  positions  deter- 
mined by  absolute  positioning  methods  may  be  used  as  mapping  control  points,  datum 
transformation  ties  or  in  other  mapping  related  applications.  Once  the  absolute  position 
of  a  station  is  determined  it  may  be  subsequently  used  as  a  known  station  in  a  relative 
positioning  survey  to  establish  the  positions  of  additional  stations. 

The  tool  used  to  produce  the  absolute  position  solutions  was  the  Geodetic  Ab- 
solute Sequential  Positioning  (GASP)  System  of  programs  developed  by  the  Defense 
Mapping  Agency  (DMA)  [Ref.  1:  p.  1,  Ref.  2:  p.  487].  GASP  estimates  the  absolute 
position  of  a  point  in  the  geocentric,  earth-fixed,  cartesian,  Conventional  Terrestrial 
(CT)  coordinate  system  in  the  World  Geodetic  System  1984  (WGS  84)  reference  frame 
(See  Figure  1).  One  axis  of  the  CT  system  passes  through  the  intersection  of  the 
Greenwich  meridian  and  the  equatorial  plane.  The  third  axis  passes  through  the  Con- 
ventional International  Origin  (CIO)  which  is  the  average  position  of  the  earth's  rota- 
tional pole  for  the  years  1900  to  1905.  The  second  axis  is  orthogonal  to  the  first  and 
third  axes  in  a  right-hand  sense.  [Ref.  3:  Sec.5.6] 
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Figure  1.      The  Conventional  Terrestrial  Coordinate  System 


DMA  has  been  processing  GPS  data  collected  with  the  TI  4100  with  GASP  for 
the  last  few  years  and  it  has  demonstrated  that  geodetic  quality  point  positions  are  rou- 
tinely achievable.  This  means  that  a  position  solution  with  a  standard  deviation  on  each 
component  of  less  than  one  meter  can  be  expected,  usually  with  data  collected  over  a 
span  of  about  four  hours  at  a  data  collection  rate  of  30  seconds.  Because  accuracy  re- 
sults using  the  Tl  4100  data  with  GASP  are  well  documented,  a  benchmark  had  already 
been  established  that  would  provide  a  reference  for  the  results  we  would  obtain.  [Ref. 
1:  p.  5] 

a.     Thesis  Objectives 

Before  this  project,  GASP  had  the  ability  to  process  data  collected  with  the 
TI  4100  receiver,  using  the  Floating  point,  Integer,  Character,  ASCII  (or  FICA)  format 
as  the  input  format.  One  of  the  goals  of  this  thesis  was  to  modify  GASP  such  that  the 
programs  would  accept  the  measurements  of  any  GPS  receiver  type.  The  Receiver 
INdependent  EXchange  (RINEX)  format  was  chosen  as  a  second  data  format  through 
which  receiver  measurements  could  be  entered  into  the  program  [Ref.  4].  The  RINEX 
format  was  selected  because  it  is  one  of  the  more  widely  recognized  and  accepted  of  the 
exchange  formats  currently  in  use  and  because  both  the  Trimble  and  Ashtech  data 
processing  packages  have  programs  that  will  convert  their  data  to  the  RINEX  format. 
It  was  also  desired  that  the  GASP  programs  retain  the  capability  to  process  data  in  the 
FICA  format  for  which  GASP  was  originally  written. 

In  addition  to  modifying  the  GASP  program  to  accept  the  data  from  the 
Trimble  and  Ashtech  receivers  other  objectives  of  this  thesis  are: 

•  To  compare  GPS  position  solutions  produced  for  the  TI  4100,  Trimble  4000ST, 
and  Ashtech  LD  XII  GPS  receivers  to  the  solutions  obtained  from  an  independent 
method.l 

•  To  compare  the  receiver  position  solutions  for  the  broadcast  ephemeris  versus  the 
precise  ephemeris 

•  To  examine  repeatability  (i.e.,  consistency)  of  results 

•  To  examine  the  effects  of  Selective  Availability  on  the  position  solution  when  using 
the  broadcast  ephemeris 


1  The  Transit  Doppler  positioning  system  was  the  independent  method  used  in  this  compar- 
ison and  the  positions  derived  from  this  method  are  taken  as  the  "true*  positions. 


B.     BACKGROUND 

1.    The  Global  Positioning  System  (GPS)  -  Fundamental  Concepts 

There  are  a  number  of  satellite  based  positioning  systems  being  used  to  establish 
the  position  of  an  observer  on  or  near  the  surface  of  the  earth.  The  Global  Positioning 
System  (GPS)  is  one  such  system  [Ref.  3:  Sec.  3.0].  The  NAVSTAR  Global  Positioning 
System  is  a  passive  navigation  and  satellite  positioning  system  operated  by  the  Depart- 
ment of  Defense  (DoD).  It  was  developed  to  afford  the  user  instantaneous  three- 
dimensional  position  information  anywhere  in  the  world.  GPS  exploits  simultaneously 
received  radio  frequency  signals  to  determine  range  measurements  between  satellites  and 
earth  based  receivers.  These  measurements  along  with  a  knowledge  of  the  satellite  posi- 
tions can  be  used  to  solve  for  the  receiver  coordinates. 

Through  recent  advances  and  upgrades  in  hardware,  software  and  the  develop- 
ment of  new  techniques  in  data  processing,  the  systems  applications  have  expanded 
dramatically.  Particularly  in  the  areas  of  high  precision  surveying  and  crustal  deforma- 
tion studies,  GPS  has  rapidly  become  competitive  with  other  positioning  systems  and 
techniques  [Ref.  5]. 

The  strength  of  the  GPS  system  in  geodetic  work  lies  in  its  relative  affordability, 
portabilty,  ease  of  operation,  and  high  accuracy  in  comparison  to  other  high  precision 
positioning  systems  [Ref.  6].  For  example,  the  Transit  Doppler  satellite  based  position- 
ing system  requires  two  to  four  days  of  data  collection  to  produce  a  geodetic  quality 
point  position.  GPS,  on  the  other  hand,  is  able  to  produce  a  geodetic  quality  point  po- 
sition in  only  four  hours  of  data  collection  [Ref.  1:  p.  5].  From  a  mapping  and  charting 
perspective,  the  fact  that  it  provides  position  determinations  in  a  unified  coordinate 
system  may  be  its  most  important  feature. 
a.     GPS  System  Components 

GPS  consists  of  three  primary  segments:  the  satellites,  the  ground  control, 
and  the  users. 

(1)  The  Satellite  Segment.  The  GPS  satellite  constellation,  as  of  mid 
1991,  consists  of  16  satellites  deployed  in  high  earth  orbit  (about  20000  km  altitude) 
configured  such  that  a  minimum  of  four  satellites  are  visible  to  the  user  at  a  given  time 
with  from  five  to  seven  satellites  typically  available.  The  satellites  are  arranged  in  six 
orbital  planes  and  have  an  orbital  period  of  about  12  hours.  It  is  proposed  that  24  sat- 
ellites eventually  be  deployed  to  ensure  that  a  rninimum  of  six  satellites  will  be  visible  to 
the  observer. 


Each  satellite  transmits  two  radio  frequency  signals  in  the  L-band, 
one  at  a  frequency  of  1575  MHz  (termed  LI)  the  other  at  1227  MHz  (L2).  The  LI  fre- 
quency is  modulated  by  the  Coarse/Acquisition  (C/A)  code,  the  Precise  (P)  code  and  the 
Navigation  Message.  The  L2  frequency  is  modulated  only  by  the  P-code  and  the  Navi- 
gation Message. 

The  C/A  and  P-codes  are  binary,  pseudo-random  noise  (PRN)  codes. 
Pseudo-random  noise  codes  resemble  true  random  noise  except  that  random  noise  car- 
ries no  information  whereas  PRN  codes  are  generated  to  be  predictable  and  carry  in- 
formation. In  this  case,  the  information  is  used  to  determine  the  time  of  signal  transit 
from  satellite  to  receiver.  Identical  codes  are  generated  by  both  the  satellites  and  the 
receiver  and  cross-correlations  (determination  of  the  scalar  product  of  the  code  sequence 
with  a  time  delayed  copy  of  itself)  of  the  incoming  satellite  codes  with  the  receiver  gen- 
erated replicas  are  performed  to  determine  the  time  reading  of  the  transmitter  clock. 

The  C/A  and  P-code  modulations  shift  the  phase  of  the  carrier  L- 
band  frequencies  by  180  degrees.  They  are  essentially  a  sequence  of  positive  and  negative 
ones  superimposed  onto  the  carriers  at  frequencies  of  1.023  MHz  and  10.23  MHz  re- 
spectively. If  the  code  value  is  minus  one  the  carrier  phase  is  shifted,  if  plus  one  there 
is  no  effect  on  the  carrier  signal  (See  Figure  2). 

The  C/A-code  sequence  repeats  every  millisecond  yielding  a  C/A  code 
cycle  of  300  kilometers.  There  are  210  C/A  code  chips  in  one  code  cycle  so  that  the  length 
of  one  code  chip  corresponds  to  about  300  meters.  Since  it  is  possible  for  receivers  to 
measure  fractions  of  chips,  the  C/A  code  may  be  used  as  a  medium  accuracy  navigation 
signal.  C/A  codes  are  exclusive  to  a  particular  satellite  making  it  possible  to  distinguish 
between  signals  received  simultaneously  from  the  satellites. 

The  P-code  sequence  repeats  every  267  days  and  is  subdivided  into 
38  seven  day  segments.  Each  satellite  is  assigned  a  one  week  segment  of  the  code.  Thus 
all  satellites  can  transmit  on  the  same  frequency  and  still  be  distinguished  from  each 
other.  This  weekly  subdivision  creates  an  identification  system  based  on  the  PRN  seg- 
ment assigned  to  a  particular  satellite.  If,  for  example,  a  satellite  is  assigned  the  seventh 
weekly  segment  of  the  code  sequence,  it  is  identified  as  PRN  7.  Codes  are  initialized 
once  per  week  at  Saturday  midnight.  Since  the  length  of  a  P-code  chip  corresponds  to 
about  30  meters,  the  P-code  supplies  a  more  precise  measurement  than  the  C/A-code. 

Critical  satellite  hardware  components  are  the  on-board  atomic 
oscillators  (highly  stable  and  precise  cesium  and  rubidium  clocks)  which  control  the 
generation  of  the  carrier  frequencies  and  code  modulations.  The  signals  are  coherently 
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Figure  2.      Code  Modulation 


generated  from  the  same  fundamental  frequency,  designated  f0  —  10.23  MHz.  The  two 
carrier  frequencies  are  multiples  of  the  fundamental  with,/^,  =  154^  andfu  =  120  f0  The 
C/A-code  modulation  has  a  chip  rate  one-tenth  the  fundamental  frequency  (fj  10)  and 
the  P-code  a  chip  rate  equal  to  the  fundamental  frequency  (f0). 

The  Navigation  Message  is  a  low  data  rate  message  (broadcast  at  50 
bits  per  second)  that  contains  information  on  satellite  health,  satellite  clocks  and  the 
Broadcast  Ephemerides.  The  contents  of  the  navigation  message  are  updated  hourly. 

Before  the  receiver  position  can  be  determined,  the  positions  of  the 
orbiting  satellites  must  be  known.  The  broadcast  ephemerides  are  a  set  of  predicted  or 
extrapolated  orbital  parameters  which  define  the  satellites  position  with  time  and  are 
used  principally  in  real-time  positioning.   Satellite  tracking  information,  received  by  five 


ground  tracking  stations  (identified  in  the  next  section),  is  used  in  a  least  squares  ad- 
justment to  estimate  the  satellite  positions.  The  positions  of  the  ground  tracking  stations 
are  held  fixed  in  the  estimation  process.  The  broadcast  ephemeris  information  is  valid 
only  over  a  specified  time  interval  (about  six  hours),  not  over  the  entire  orbit.  The  main 
reason  the  navigation  message  is  updated  hourly  is  so  the  broadcast  ephemeris  remains 
current. 

The  Keplerian  orbit  description  contained  in  the  broadcast  ephemeris 
is  the  means  by  which  the  conventional  terrestrial  coordinates  of  the  satellites  are  com- 
puted. Six  Keplerian  elements  and  a  reference  time  are  needed  to  completely  describe  the 
satellites  orbit  in  the  CT  coordinate  system.  Five  of  these  elements  define  the  satellite 
orbit  and  the  other  describes  the  position  of  the  satellite  in  the  orbit  as  a  function  of 
time.  The  other  broadcast  ephemeris  parameters  transmitted  in  the  Navigation  Message 
describe  the  deviations  or  perturbations  of  the  satellite  motion  from  the  smooth  ellipse 
defined  by  the  six  Keplerian  elements  refered  to  above.  For  the  complete  list  of  the 
broadcast  parameters  and  the  computations  required  to  convert  the  Keplerian  orbital 
parameters  to  coordinates  in  the  conventional  terrestrial  system,  see  Appendix  A. 

If  more  accurately  determined  satellite  orbits  are  required  than  are 
provided  in  the  broadcast  ephemeris,  a  post-computed  or  precise  ephemeris  may  be  ob- 
tained. The  Defense  Mapping  Agency  and  the  U.S.  National  Geodetic  Service  are  the 
agencies  responsible  for  generating  and  distributing  this  information  to  the  user  upon 
request.  Because  the  precise  ephemeris  is  computed  after  satellite  observations  are  made 
and  is  an  interpolation  of  satellite  position  based  on  a  least  squares  adjustment  using  ten 
ground  stations  rather  than  five,  much  better  accuracies  for  the  satellite  position  are  at- 
tainable. The  precise  ephemerides,  as  supplied  by  either  DMA  or  NGS,  differ  from  the 
broadcast  ephemerides  in  both  data  presentation  and  in  the  frequency  with  which  satel- 
lite positions  are  estimated.  The  precise  ephemeris  furnishes  satellite  positions  already 
computed  in  the  CT  coordinate  system  and  gives  the  satellite  velocities  in  the  directions 
of  the  coordinate  axes.  The  CT  satellite  positions  and  velocities  are  determined  every 
15  minutes  in  contrast  to  the  hourly  updates  of  the  broadcast  ephemerides. 

(2)  The  Ground  Control  Segment.  This  segment  consists  of  five  globally 
distributed  monitoring  stations  that  track  the  satellites  and  transmit  tracking  informa- 
tion to  the  master  control  station  in  Colorado  Springs,  Colorado.  It  is  here  that  com- 
putations are  performed  to  provide  the  broadcast  ephemeris  contained  in  the  updated 
navigation  message.  The  updated  navigation  message  is  then  uploaded  to  the  satellites 
for  broadcast  to  the  user.  The  navigation  message,  as  previously  stated,  contains  infor- 


mation  on  the  health  status  of  the  satellites,  and  information  on  the  satellites  orbits  and 
atomic  clocks.  The  five  monitoring  stations  are  located  at  Diego  Garcia,  Ascension  Is- 
land, Kwajalein,  Hawaii,  and  Colorado  Springs  (See  Figure  3). 


Figure  3.      GPS  Control  Stations 


(3)  The  User  Segment.  The  user  segment  is  composed  of  civilian  and 
military  users  exploiting  some  type  of  GPS  receiver  and  receiving  antenna  system.  GPS 
receivers  have  the  ability  to  receive  the  LI  and  L2  signals  from  a  number  of  satellites 
simultaneously  by  devoting  specific  channels  to  each  satellite  signal.  Once  the  signals  are 
channeled  they  can  be  individually  processed.  The  three  most  commonly  used  types  of 
GPS  receiver  are  the  code-correlating  receivers,  squaring  channel  receivers,  and  receivers 
that  incorporate  elements  of  both  of  these  techniques.  The  code-correlating  technique 
is  normally  employed  to  decode  the  broadcast  message  and  also  to  provide  pseudorange 


measurements  from  either  the  C/A  or  P-codes.  The  TI  4100  provides  pseudorange 
measurements  from  the  P-code  while  the  Trimble  and  Ashtech  pseudoranges  are  derived 
from  the  C/A-code.  Once  the  carrier  signal  has  been  demodulated  (i.e.  the  navigation 
message  and  PRN  codes  removed),  the  carrier  frequencies  may  be  processed  to  provide 
carrier  phase  measurements.  The  squaring  channel  technique  provides  carrier  phase 
measurements  but  cannot  be  used  to  extract  the  broadcast  message  (satellite  position 
and  clock  information  must  be  supplied  externally)  nor  provide  pseudoranges  (for  details 
on  the  pseudorange  and  carrier  phase  measurements  see  the  next  section). 

Essentially,  code-correlating  receivers  operate  by  internally  generating 
a  replica  of  the  incoming  code  modulated  satellite  signal.  The  start  of  the  code  replica 
will  be  offset  due  to  the  propagation  delay  between  satellite  and  receiver.  The  replica  is 
then  incrementally  shifted  or  cross-correlated  with  the  incoming  coded  signal  and  the 
clock  controlling  the  generation  of  the  replica  signal  is  corrected  to  reflect  the  shifts. 
Once  the  replica  is  aligned  to  the  incoming  signal,  it  stays  locked  to  it.  At  this  point,  the 
replica  code  generator  clock  reads  identically  to  the  satellite  clock  and  the  signal  trans- 
mission time  is  determined  by  differencing  this  time  from  the  time  recorded  on  a  receiver 
clock  on  GPS  time.  This  time  difference  multiplied  by  the  speed  of  signal  propagation 
gives  the  pseudorange.  Once  code  lock  is  attained,  the  code  may  be  removed  from  the 
incoming  signal  and  the  satellite  navigation  message  extracted.  The  signal,  now  demod- 
ulated, may  at  this  point  be  used  for  processing  carrier  phase  measurements. 

Squaring  channel  receivers  operate  by  squaring  the  incoming  satellite 
signal  effectively  removing  any  phase  reversals  resulting  from  code  modulations  im- 
pressed on  the  signal  at  the  satellite.  This  produces  a  signal  that  is  double  the  frequency 
of  the  incoming  signal  (See  Figure  4).  The  phase  difference  between  this  squared  signal 
and  an  internally  generated  replica  is  still  easily  determined  however.  The  accumulated 
phase  difference  (whole  and  fractional  cycles)  is  continuously  counted  and  from  this  the 
carrier  based  change  in  range  is  determined  by  multiplying  by  the  wavelength  of  the 
squared  signal. 
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Figure  4.      Squaring  Channel 

GPS  receivers  are  equipped  with  quartz  crystal  oscillators  rather  than 
the  more  precise,  stable,  and  expensive  atomic  oscillators  used  aboard  the  satellites. 
These  oscillators  do  however  provide  more  than  sufficient  precision  and  short  term  sta- 
bility for  the  typical  length  of  a  surveying  session.  Like  the  satellite  clocks,  the  receiver 
clocks  control  signal  generation. 

Internal  receiver  software  estimates  the  receiver  position  by  employ- 
ing least  squares  or  Kalman  filtering  algorithms  using  the  broadcast  ephemeris  satellite 
positions  and  pseudorange  or  carrier  phase  measurements.  All  receivers  estimate  the 
receiver  clock  bias  in  addition  to  the  postion  parameters  and  then  use  this  value  to  reset 
the  receiver  clock  to  GPS  time. 
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b.     GPS  Operating  Principles  -  Emphasis  on  Point  Positioning 

(J)  GPS  Observables  and  Observation  Equations.  There  are  two  types  of 
observation  that  may  be  exploited  to  determine  positions,  the  pseudorange  measurement 
and  the  carrier  phase  measurement.  For  point  positioning  particularly  with  GASP,  the 
carrier  phase  is  of  primary  interest  and  so  is  emphasized  here.  Pseudorange  measure- 
ments to  the  extent  they  are  used  by  GASP,  serve  only  to  monitor  the  validity  of  the 
carrier  phase  measurements  and  are  not  used  in  the  actual  position  estimation.  The 
fundamentals  of  pseudorange  positioning  are  briefly  touched  on,  however,  in  an  effort 
to  convey  some  general  concepts  and  introduce  some  common  terms. 

(2)  Pseudorange.  Under  ideal  conditions  (i.e.,  precisely  known  satellite 
location,  signal  propagation  speed,  and  perfectly  synchronized  satellite  and  receiver 
clocks)  the  travel  time  of  the  signal  would  be  given  by 

where  tr  is  the  time  of  signal  reception  at  the  receiver  and  tx  is  the  time  of  satellite  signal 
transmission.  This  time  difference  would  then  be  converted  to  distance  units  (p)  by 
multiplying  by  the  appropriate  speed  of  signal  propagation  (c),  where  p  —  c-c.  This 
range  measurement,  along  with  the  range  measurements  of  two  additional  satellites  and 
the  known  locations  of  the  three  satellites,  could  then  produce  a  system  of  equations 
that  can  be  solved  to  uniquely  determine  the  receiver  positions. 

In  the  real  world,  we  must  be  concerned  with  error  sources  and  their 
effects.  A  major  source  of  error  in  any  GPS  measurement  is  the  inability  of  the  receiver 
and  satellite  clocks  to  maintain  alignment  with  the  reference  time  standard  (GPS  time). 
Of  these  two  types  of  clock  error,  that  associated  with  the  receiver  dominates  due  to  the 
lower  precision  and  stability  of  these  clocks.  The  Allan  variance  is  often  cited  as  an  in- 
dication of  oscillator  stability.  The  Allan  variance  values  for  receiver  crystal  quartz 
oscillators  are  typically  several  orders  of  magnitude  larger  than  the  satellite  atomic 
oscillators.  Satellite  clocks  take  longer  to  warm-up  and  to  stabilize  initially  but  they 
maintain  stability  much  longer  than  do  receiver  clocks. 

To  address  the  problem  of  clock  errors,  we  introduce  a  combined 
satellite  and  receiver,  first  order  clock  correction  term  as  a  fourth  parameter  in  our  sys- 
tem of  equations.  Expressing  the  receiver  clock  offset  from  GPS  time  as 

tr  =  tr-t 
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where, 

tr  is  the  receiver  clock,  time 
t  is  GPS  time 

and  the  satellite  clock  offset  from  GPS  time  as 


Tx  =  *x  ~  * 


where, 

tx  is  the  satellite  clock  time 

the  combined  clock  correction  term  becomes 

dx  =  rr  -  rx 

Because  of  clock  and  other  types  of  errors,  the  measured  range  between  satellite  and 
receiver  is  not  the  true  range  but  is  rather  a  biased  range  referred  to  as  the  pseudorange. 
Designating  the  pseudorange  as  p,  the  system  of  observation  equations  is  represented 
by: 

Pr,  =  Pt  +  c  dTi 

with  i  greater  than  or  equal  to  4, 

where  p„  the  true  range  from  the  ith  satellite,  is  given  by 

Pi  =  IK  -  xrf  +  (ySi  -yrf  +  (zSi  -  zr)2],/2  ~\XSl-XT\ 

In  this  equation, 

x.  is  the  x  coordinate  of  the  ith  satellite 
'i 

y,t  is  the  y  coordinate  of  the  ith  satellite 
2..  is  the  z  coordinate  of  the  ith  satellite  and 
x„y„  zr  are  the  unknown  receiver  coordinates. 

The  addition  of  this  fourth  parameter  requires  a  minimum  of  four  observation  equations 
to  solve  for  the  receiver  coordinates  and  the  clock  correction.  If  xp  is  now  introduced  to 
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express  the  time  delay  associated  with  other  error  sources  the  pseudorange  equation 
becomes: 

Prt  =  Pi  +  c  {dxt  +  rp) 

The  pseudorange  is  the  observable  most  often  used  to  obtain  a  navi- 
gation solution  where  submeter  level  accuracy  is  not  required.  Three  dimensional  navi- 
gation solutions  currently  offer  best  case  positional  accuracy  at  16  meters.  While 
certainly  suitable  for  navigating  a  vessel,  a  position  in  error  by  16  meters  could  not  be 
used  as  a  reference  for  establishing  or  extending  mapping  control.  For  high  precision 
relative  or  point  positioning,  phase  observables  are  needed.  Because  the  carrier  phase 
wavelength  is  shorter  than  the  code  modulations  to  the  carrier  signal  (20  centimeters  for 
LI  as  opposed  to  30  meters  for  the  P-code),  the  carrier  signal  may  be  used  to  provide  a 
more  precise  distance  measurement. 

(3)  Carrier  Phase.  It  is  possible  to  obtain  distance  (or  more  exactly 
change  of  distance)  information  by  measuring  the  phase  of  the  carrier  signal.  The  dif- 
ference between  the  phase  of  a  receiver  generated  carrier  signal  and  the  incoming, 
Doppler  shifted,  satellite  carrier  signal  gives  the  carrier  beat  phase  observable 

/ 

<Pmeas  =  ~T  P  +fdr  ~fXP 

where,  <f>  is  in  cycles  and  /is  the  signal  frequency. 

The  total  continuous  carrier  beat  phase  measurement  consists  of  an 
accumulating  count  of  the  whole  and  fractional  difference  in  cycles  since  the  time  of 
signal  acquisition.  What  is  not  recorded  is  an  unknown  number  of  integer  cycles  at  initial 
signal  reception  (the  whole  number  of  cycles  between  the  receiver  and  satellite).  The 
unknown  number  of  cycles  is  referred  to  as  the  cycle  or  integer  ambiguity.  The  total 
continuous  carrier  phase  can  be  expressed 

4>  total  =  <t>meas  +  N(h) 

where  N(t0)  is  the  integer  ambiguity.     A  carrier  phase  equation  comparable  to  the 
pseudorange  equation  can  be  written, 


meas, 


P<t>,  =  *■<!>* 

=  pl-c{dT,  +  Xp)  +  XN(t0)l 
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where,  p^  is  the  phase  biased  pseudorange, 

tp  is  the  signal  time  delay  associated  with  other  error  sources  and 
X  —  elf  is  the  signal  wavelength. 


See  Figure  5  for 

a  graphic  representation  of  the  carrier  phase  concept. 
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Figure  5.      Components  of  Carrier  Phase  Measurements 

c.     Measurement  Errors  and  Error  Models 

There  are  two  types  of  errors  that  effect  GPS  measurements  and  degrade  the 
accuracy  of  the  estimated  position,  systematic  errors  and  random  errors.  Systematic 
errors  are  errors  resulting  from  a  predictable  source.  These  errors  are  typically  of  con- 
stant magnitude  under  a  set  of  given  conditions.  When  the  cause  is  understood,  these 
errors  can  be  removed  by  a.  correction  model  or  a  change  in  observation  procedure. 
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Measurement  biases  resulting  from  systematic  errors  produce  offsets  between  the  values 
that  are  observed  and  the  "true"  value.  After  systematic  errors  are  eliminated  from 
measurements,  random  errors  remain.  These  errors  are  defined  by  repeated  measure- 
ments. A  repeated  measurement  will  not  yield  identical  values  each  time  due  to  the  oc- 
currence of  random  errors.  They  indicate  how  well  a  measurement  may  be  repeated  and 
characterize  the  precision  of  a  measurement.  If  they  are  the  only  error  source,  they  are 
a  measure  of  accuracy.  Random  errors  are  usually  small  and  the  probability  of  a  positive 
or  negative  error  of  a  given  magnitude  are  the  same.  Random  errors  are  dealt  with  in  a 
least  squares  adjustment.  [Ref.  7,  8] 

GPS  systematic  errors  can  be  divided  into  three  general  categories:  satellite 
errors,  station  errors,  and  observation  dependent  errors.  Satellite  errors  include  errors  in 
the  computed  satellite  ephemeris  or  in  the  model  for  the  satellite  clock  provided  in  the 
navigation  message.  Station  errors  occur  as  a  result  of  errors  in  the  receiver  clock.  Ob- 
servation dependent  errors  are  related  to  factors  that  influence  the  speed  of  signal 
propagation,  ambiguities  in  phase  observables,  etc. 

(I)  Satellite  and  Receiver  Clock  Errors.  Inability  to  perfectly  synchro- 
nize the  receiver  clock  and  the  more  tightly  controlled  satellite  clock  to  a  standard  ref- 
erence (or  GPS)  time  produces  time  offsets.  If,  for  instance,  there  exists  a  one 
microsecond  timing  misalignment  between  satellite  time  and  receiver  time,  a  300  meter 
range  bias  will  result.  Additionally,  the  clocks  may  be  misaligned  in  frequency  (a- fre- 
quency bias)  or  the  frequencies  may  change  over  time  (frequency  drift).  Any  of  these 
will  produce  clock  errors.  Random  error  (or  noise)  also  effects  the  time  measurements. 
The  timing  behavior  of  the  satellite  clocks  is  carefully  monitored  by 
the  ground  stations  and  the  drift  of  the  clocks  from  standard  GPS  time  determined.  The 
amount  of  the  satellite  time  offset,  frequency  offset,  and  frequency  drift  are  then  ex- 
pressed as  coefficients  of  a  second-order  polynomial  that  are  transmitted  in  the  naviga- 
tion message.  The  polynomial  coefficients  are  determined  via  a  least  squares  adjustment 
performed  at  the  master  control  station  in  Colorado  Springs.  The  polynomial  has  the 
form, 

T*  -  %  +  «i  ('  ~  'o)  +  ai  ('  ~  'o)2 

where,  /„  is  some  reference  epoch, 

4,  is  the  satellite  clock  time  offset, 
ax  is  the  frequency  offset,  and 
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a2  is  the  frequency  drift.    In  this  manner,  satellite  clock  synchronization  to  GPS 
time  is  typically  maintained  to  within  20  nanoseconds. 

GPS  receiver  clock  errors  must  be  modeled  if  the  user  desires  a  high 
accuracy  solution.  Receiver  clock  errors  may  be  treated  in  a  fashion  similiar  to  that  for 
the  satellite  clocks.  The  polynomial  coefficients  may  be  estimated  as  additional  param- 
eters in  a  least  squares  adjustment  along  with  the  receiver  coordinates. 

(2)  Orbit  Errors.  The  positions  of  the  satellites  with  time  are  well  known 
but,  due  to  forces  acting  on  the  satellites  that  may  not  be  adequately  modeled,  are  not 
perfectly  determined.  The  error  in  satellite  position  propagates  to  contribute  to  error  in 
the  receiver  position.  The  magnitude  of  this  error  depends  on  whether  the  broadcast 
(predicted)  ephemeris  or  the  precise  (post-fit)  ephemeris  is  used  to  determine  the  satellite 
positions. 

Improving  the  models  used  to  describe  the  forces  acting  on  the  satel- 
lite is  one  approach  to  obtaining  more  accurate  satellite  positions.  Another  is  to  include 
parametric  models  for  the  forces  as  part  of  the  orbit  estimation  process  performed  by  the 
control  segment.  Depending  on  the  application,  the  effects  of  orbit  errors  may  be  ig- 
nored altogether  or  the  data  may  be  differenced  to  reduce  or  eliminate  these  effects. 

(3)  Observation  Dependent  Errors.  The  speed  of  signal  propagation  is 
influenced  by  many  factors  including;  signal  interaction  with  the  ionosphere  and 
troposphere,  relativistic  effects  on  the  signal,  etc. 

(a)  Signal  Interaction  with  the  Ionosphere—  Signal  interaction 
with  the  free  electrons  found  in  the  ionosphere  produces  a  change  in  path  length  that 
may  be  on  the  order  of  tens  of  meters.  The  free  electrons  are  released  from  gas  molecules 
ionized  by  incoming  solar  ultra-violet  radiation.  Any  condition  that  acts  to  release  more 
electrons,  such  as  increasing  the  amount  of  incoming  solar  radiation  (midday  or  during 
increased  sunspot  activity),  will  correspondingly  lengthen  the  signal  path.  Satellite- 
receiver  geometry  also  plays  a  role  as  the  total  number  of  free  electrons  along  the  path 
is  a  function  of  the  distance  that  the  signal  travels  through  the  ionosphere.  Thus  when 
a  satellite  is  near  the  horizon  the  signal  encounters  more  electrons  than  when  near  the 
zenith.  This  ionospheric  effect  is  frequency  dependent  and  is  inversely  proportional  to 
the  square  of  the  signal  frequency.  A  comparison  of  measurements  on  LI  and  L2  may 
be  used  to  derive  a  dual  frequency  correction  given  by 

dPL\ F 
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where  T  is  defined  by 

flA     \2 


r  = 


In  J 


For  single  frequency  receivers,  an  ionospheric  model  must  be  used  in  place  of  the  dual 
frequency  correction  or  ignored  altogether.  [Ref.  2:  p.  489] 

(b)  Tropospheric  Interaction—  Tropospheric  effects  are  a  result 
of  refraction  in  the  neutral  atmosphere  and  are  not  frequency  dependent.  This  effect  may 
be  separated  into  two  components:  the  dry  component,  and  the  wet  component.  The  dry 
component  comprises  about  90%  of  the  total  effect  and  is  a  function  of  the  surface  at- 
mospheric pressure.  It  is  approximated  in  the  vertical  by 

DTC=  2.27  x  \QT\mlmb)  P0 

where  P0  (the  atmospheric  pressure)  is  in  millibars  and  DTC  (the  dry  term  range  con- 
tribution) is  in  meters.  In  the  zenith  direction,  this  corresponds  to  a  maximum  range 
bias  of  about  2.5  meters.  The  DTC  increases  with  decreasing  satellite  elevation  angle  and 
at  five  degrees  above  the  horizon  it  ranges  from  20-30  meters. 

Estimating  the  effects  of  the  wet  component  is  a  more  compli- 
cated proposition.  It  depends  on  the  total  water  vapor  content  along  the  signal  path  and 
hence  on  the  temperature,  pressure,  and  humidity.  Surface  temperature,  pressure  and 
humidity  values  are  used  to  estimate  the  magnitude  of  the  integrated  effect  by  employing 
an  atmospheric  model  such  as  the  Hopfield  or  Chao  models  [Ref.  2:  p.  490].  If  obser- 
vations from  water  vapor  radiometers  are  available,  they  will  most  accurately  profile 
atmospheric  conditions  along  the  path. 

(c)  Integer  Ambiguity—  In  differential  positioning  applications, 
missing  whole  cycles  (refered  to  as  integer  ambiguities)  in  the  carrier  phase  measure- 
ments must  be  resolved  if  we  are  to  fully  exploit  the  more  precise  nature  of  phase 
measurements.  Efforts  must  be  made  to  account  for  any  whole  cycles  not  recorded  in 
the  observations.  This  may  be  the  result  of  a  loss  of  phase  lock  between  receiver  and 
satellite  due  to  some  obstruction  to  the  signal.  Also  the  initial  number  of  whole  cycles 
between  satellite  and  receiver  (the  integer  ambiguity)  must  be  determined.  Many  tech- 
niques have  been  developed  to  resolve  the  absense  of  any  cycles  from  the  measurements. 
In  point  positioning  with  GASP,  resolution  of  the  integer  ambiguity  is  not  a  problem 
as  we  shall  see. 
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Other  effects  that  must  be  accounted  for  include:  relativistic 
effects,  earth  rotation  (i.e.,  receiver  position  in  motion)  during  time  of  signal  trans- 
mission, and  offset  of  the  satellite  transmitting  antenna  position  from  the  position  of  the 
satellite  center  of  mass  (the  position  given  by  the  ephemeris).  Corrections  to  account  for 
the  effects  of  these  error  sources  are  applied  to  the  phase  biased  pseudorange  and  the 
pseudorange  observables  before  the  GASP  model  is  formed  and  the  estimation  algorithm 
is  executed.  Details  on  the  models  employed  to  compute  these  corrections  are  given  in 
the  chapter  on  GASP  processing. 

(4)  Comments  on  Selective  Availability.  Selective  Availability  (SA)  is  the 
intentional  degradation  of  the  position  solution  available  to  a  select  segment  of  the  user 
community  (i.e.,  most  civilians  and  unfriendly  military)  operating  in  a  real-time  naviga- 
tion mode.  It  may  be  implemented  by  the  DoD  at  their  discretion.  The  solution  degra- 
dation is  accomplished  by  broadcasting  inaccurate  positions  for  the  orbiting  satellites 
or  by  dithering  the  satellite  clock  so  that  an  inaccurate  signal  transmission  time  is  ob- 
tained. This  will  not  greatly  effect  relative  positioning  applications  since  receiver  differ- 
encing schemes  remove  these  errors.  For  point  positioning  applications  that  rely  on  the 
post-computed  precise  ephemeris,  only  the  clock  dither  will  effect  the  solutions. 
d.     Geometric  Effects 

In  addition  to  the  errors  just  presented,  the  configuration  of  the  satellites 
relative  to  the  receiver  also  has  an  effect  on  the  determination  of  the  receiver  position. 
This  geometric  influence  is  referred  to  as  the  Dilution  of  Precision  or  DOP  factor.  It  is 
roughly  represented  by  the  ratio  of  the  positioning  accuracy  (ap)  to  the  measurement 
accuracy  (a0)  or 

DOPx-rf- 

ao 

Actually  there  are  a  number  of  DOP  factors.  The  most  commonly  referred  to  of  these, 
the  Geometric  Dilution  of  Precision  (GDOP),  is  more  precisely  defined  as  the  square 
root  of  the  trace  of  the  covariance  matrix  (the  trace  being  the  sum  of  the  four  diagonal 
elements;  the  three  position  variances  and  the  time  variance).  The  DOP  is  a  measure 
of  the  geometric  strength  of  the  satellite  configuration  that  changes  with  time.  The  time 
of  most  favorable  observation  for  obtaining  a  user  position  is  when  the  DOP  is  small 
(with  a  DOP  of  less  than  five  preferred)  and  falling.  Since  the  number  and  the  positions 
of  the  satellites  visible  to  the  observer  changes  over  time,  many  possible  geometric  sce- 
narios are  available  over  a  tracking  session.   A  GPS  survey  should  be  planned  to  take 
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best  advantage  of  changing  geometric  scenarios.     For  point  positioning,  a  geodetic 
quality  position  solution  can  usually  be  obtained  with  approximately  four  hours  of  data 
collection.  During  times  of  very  favorable  geometry,  a  geodetic  quality  point  position 
may  be  obtained  in  a  shorter  time. 
e.     Differencing  Techniques 

By  forming  linear  combinations  of  the  basic  pseudorange  or  carrier  phase 
equations,  a  number  of  error  sources  common  to  the  measurements  being  differenced 
will  cancel  or  be  greatly  reduced.  The  accuracy  of  the  computed  position  can  be  signif- 
icantly improved  by  employing  differencing  because  many  sources  of  measurement  error 
are  removed  or  reduced.  The  three  types  of  differencing  combinations  often  used  are: 
differencing  between  two  satellites,  differencing  between  two  receivers,  and  differencing 
between  two  epochs  (time  periods).  These  are  referred  to  as  single  differences.  The  dif- 
ferencing schemes  involving  the  carrier  phase  equations  are  the  only  ones  presented  here 
since  they  are  used  to  form  the  GASP  model.  Differencing  schemes  that  involve  the 
pseudorange  equation  are  ignored  since  pseudorange  observations  do  not  contribute  to 
the  position  estimation  except  as  a  means  of  screening  the  carrier  phase  observations. 
Also,  differencing  techniques  involving  multiple  receivers  are  neglected  in  this  treatment 
except  to  state  what  types  exist  and  the  what  errors  they  effect. 

Between-epoch  single  differencing  involves  differencing  two  equations  of  the 
instantaneous  carrier  phase  for  one  satellite  and  one  receiver.  The  carrier  phase 
equation  at  epoch  1  is  described  by 

P*,  =  Pi  ~  c  (<*Ti  +  TPl)  +  x  N(t0) 
The  carrier  phase  equation  at  epoch  2  is  given  by 

P<t>3  =  Pi  ~  c  (d?2  +  *fc)  +  ^  ^('o) 
Differencing  then  yields 

P*!  -P<t>2  =  (Pi~  Pi)  ~  c  (dri  -  d?i  +  TPl  -  Tft) 
Introducing  3  as  the  difference  notation  the  equation  may  be  written 

bp^  =  6p-c  {ddx  -  Stp) 

See  Figure  6  for  a  graphic  depiction  of  between-epoch  single  differencing. 
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Between-Epoch   Single    Differencing 


Satellite  a 


Figure  6.      Between-Epoch  Single  Differences:     Differencing  Range  Equations  of 
One  Satellite  Over  Two  Consecutive  Measurement  Epochs 

The  advantage  of  this  scheme  is  that  the  initial  integer  ambiguity  is  removed. 

Similiar  manipulations  can  be  performed  to  yield  between-satellite  single 
differences  (See  Figure  7).  For  two  satellites  represented  by  the  superscripts  a  and  /?, 
phase  equations  are  given  by, 
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p;  =  pa-c(jza  +  Tap)  +  tN(t0y 


and, 


P* 


l-if-c{dS  +  4)  +  XN{tf 


The  between-satellite  single  difference  equation  for  the  carrier  phase  is 

P+-  P*l  =  P*  -  PP  -c(ctra  -dx"  +  rap-4)  +  HN(toy  -  N(t0f) 


or, 


dp,p  =  dp  +  c  (drx  -  dxp)  +  X  dN 


Between-Satellite    Single    Differencing 


Satellite  0 


Satellite  a. 


Figure  7.      Between-Satellite  Single  Differences:     Differencing  Range  Equations  of 
Two  Satellites  Over  One  Measurement  Epoch 
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For  between-satellite  single  differencing,  receiver  clock  errors  are  removed  or  reduced. 

It  is  also  possible  to  form  between-receiver  single  difference  equations  to 
remove  satellite  clock  and  orbit  errors.  This  is  important  in  relative  positioning  but 
cannot  be  used  in  point  positioning. 

We  can  now  form  double  difference  equations  from  the  single  difference 
equations.  For  example,  a  satellite-receiver  double  difference  may  be  formed  by  differ- 
encing two  between-satellite  single  differences  (involving  the  same  pair  of  satellites)  ever 
two  receivers.  This  scheme  removes  or  reduces  the  effects  of  errors  associated  with  sat- 
ellite and  receiver  clocks. 

Receiver-time  double  differences  and  satellite-time  double  differences  may 
be  formed  in  a  similiar  manner.  Receiver-time  double  differencing  removes  or  reduces  the 
effects  of  errors  associated  with  satellite  clocks  and  eliminates  integer  ambiguities. 
Satellite-time  double  differencing  (essentially  the  GASP  model,  see  Chap.  Ill  for  more 
details)  removes  or  reduces  the  effects  of  errors  associated  the  receiver  clock  and  elimi- 
nates integer  ambiguities. 

It  is  also  possible  to  construct  an  equation  for  a  receiver-satellite-time  triple 
difference  observable.  This  gives  the  change  in  the  receiver-satellite  double  difference 
from  one  epoch  to  the  following  epoch.  With  triple  differencing,  in  addition  to  the  can- 
cellation of  integer  ambiguities  for  carrier  phase  measurements,  all  clock  and  satellite 
orbit  errors  are  removed.  Triple  differencing  is  widely  used  in  relative  positioning  appli- 
cations [Ref.  9].  This  technique  cannot  be  employed  in  point  positioning. 

The  disadvantages  to  differencing  are  that  the  number  of  observations  has 
been  reduced  and  mathematical  correlations  are  introduced  as  a  product  of  the  differ- 
encing process.  This  produces  a  weaker  solution  than  that  provided  by  not  differencing. 
Correlation  matrices  should  be  computed  to  equate  the  differenced  data  to  the  undiffer- 
enced  data  but  generally  differenced  observables  are  treated  as  uncorrelated. 

Now  that  systematic  errors  have  been  reduced  or  removed  by  error  model- 
ing (i.e.,  applying  measurement  corrections)  or  by  differencing  and  the  mathematical 
model  has  been  defined  by  forming  the  undifferenced  or  differenced  observation 
equation,  the  receiver  coordinates  and  other  parameters  may  be  estimated  by  invoking 
the  method  of  least  squares.  For  a  more  complete  explanation  of  the  use  of  the  least 
squares  method  in  GPS,  see  Appendix  B. 
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2.     Positioning  with  Transit  Doppler 

The  method  of  Transit  Doppler  positioning  is  briefly  reviewed  here.  Since  this 
method  was  used  to  provide  the  independent  reference  positions  to  which  the  GPS  de- 
rived point  positions  were  compared,  some  insight  into  Transit  Doppler  positioning  is 
necessary. 

The  U.S.  Navy  Navigation  Satellite  System  or  Transit  has  been  in  continuous 
operation  since  the  mid  1960's.  Transit  can  be  regarded  as  the  forerunner  of  GPS  in 
many  respects.  Many  of  the  ideas  and  techniques  developed  during  the  era  of  Transit 
positioning  have  been  refined  and  employed  in  GPS.  It  is  anticipated  that  Transit  will 
be  replaced  by  GPS  in  the  near  future.  Like  GPS,  Transit  consists  of  three  segments:  the 
satellites,  the  ground  control,  and  the  users. 

Six  active  Transit  satellites  are  deployed  in  circular  polar  orbit  at  an  altitude  of 
approximately  1100  kilometers.  This  height  was  selected  because  excessive  orbital  height 
would  provide  too  low  a  rate  of  change  of  Doppler  frequency.  If  the  orbit  were  too  low, 
the  Doppler  frequency  rate  of  change,  ionospheric  refraction,  and  especially  the  effects 
of  atmospheric  drag  become  too  great.  The  orbital  period  of  the  satellites  is  about  107 
minutes.  See  Figure  8  for  a  general  representation  of  the  Transit  Positioning  System. 
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Each  satellite  transmits  two  separate  frequencies  at  400  MHz  and  150  MHz. 
The  frequency  generation  is  controlled  by  a  single  highly  stable  crystal  oscillator.  The 
use  of  two  frequencies  permits  the  determination  of  an  ionospheric  refraction  correction. 
A  broadcast  message  containing  orbital  information  is  superimposed  on  the  two  carrier 
frequencies  by  phase  modulation. 

Three  major  differences  between  GPS  and  Transit  may  already  be  noted;  the 
differences  in  satellite  height,  the  use  of  only  six  satellites  as  opposed  to  24  GPS  satel- 
lites, and  the  use  of  crystal  oscillators  rather  than  the  more  precise  atomic  oscillators 
on-board  the  satellites.  By  having  only  six  satellites  available  with  fewer  than  this  actu- 
ally visible,  it  takes  much  more  time  to  acquire  an  equivalent  number  of  Transit  obser- 
vations. 

Tracking  stations  record  the  Doppler  measurements  (i.e.,  the  Doppler  shift  in 
the  frequency  transmitted  by  the  satellite)  on  each  satellite  pass.  This  information  is  re- 
layed to  the  central  processing  or  control  station  where  the  satellite  orbits  are  determined 
and  extrapolated,  then  this  ephemeris  data  is  updated  and  uploaded  to  the  satellites 
(about  every  12  hours)  for  subsequent  rebroadcast  to  the  user.  In  addition,  a  timing 
station  is  responsible  for  monitoring  the  time  signals  received  from  the  satellites  and 
adjusting  the  satellite  clocks  as  needed. 

a.     The  Basic  Principle  of  Doppler  Positioning 

The  received  frequency  will  differ  from  the  transmitted  frequency  due  the 
Doppler  effect  because  the  receiver  and  satellite  are  moving  relative  to  each  other.  If  the 
satellite  transmits  a  stable  frequency/,  then  the  frequency  at  the  receiver  is  given  by 

where, 

r  =  —  is  the  range  rate 

r  is  the  distance  or  range  between  receiver  and  satellite 

c  is  the  speed  of  signal  propagation. 

The  Doppler  frequencies  (or  Doppler  shifts)  are  measured  by  subtracting 
the  received  shifted  frequencies/  from  a  constant  receiver  reference  frequency/  (the  400 
and  150  MHz  frequencies  alluded  to  earlier).  The  time  of  closest  approach  of  the  satel- 
lite is  the  time  when/  equals/.  Integrated  Doppler  measurement  techniques  that  count 
the  number  of  accumulated  cycles  of  Doppler  shift  are  utilized  in  most  Transit  Doppler 
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receivers.  The  reason  for  this  is  that  counting  cycles  can  be  performed  more  precisely 
than  instantaneously  measuring  the  frequency. 

The  receiver  position  may  be  determined  by  continuously  counting  the 
number  of  cycles  of  the  Doppler  frequency.  Referred  to  as  the  Continuously  Integrated 
Doppler  (CID)  measurement  mode,  these  observations  along  with  accurate  positions  for 
the  satellite  can  be  used  to  establish  the  receiver  position  in  two  dimensions.  For  three- 
dimensional  positioning,  multiple  satellite  passes  must  be  observed. 

The  Doppler  counts  must  be  corrected  for  the  effects  of  atmospheric  re- 
fraction in  both  the  ionosphere  and  the  troposphere.  In  this  respect,  all  other  errors 
(relativity,  earth  rotation,  etc.)  that  effected  the  GPS  measurements  will  be  present  in  the 
Transit  Doppler  measurements  as  well.  Attempts,  such  as  error  modeling,  should  be 
made  to  minimize  their  impact. 

For  the  point  positioning  mode,  the  one  of  interest  in  this  study,  the  ob- 
servations from  multiple  satellite  passes  were  collected  with  a  single  Doppler  receiver 
over  about  a  four  day  period.  Precise  ephemerides  were  used  to  compute  the  satellite 
positions  from  which  the  CT  coordinates  of  the  site  were  determined.  The  estimated  ac- 
curacies for  a  position  solution  achievable  under  this  type  of  scenario  is  believed  to  be 
on  the  order  of  one  meter. 
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II.     EXPERIMENTAL  PROCEDURE 

A.     GENERAL 

The  Monterey  Bay  Precision  Positioning  Experiment  was  conducted  in  early  De- 
cember 1990.  Aspects  of  the  experiment  relevant  to  point  positioning  will  be  emphasized 
in  this  and  following  sections.  A  large  data  set  of  static  and  kinematic  GPS  measure- 
ments was  collected  over  four  consecutive  nights,  from  December  4  through  December 
7,  with  four  commercially  available  GPS  receivers.  Each  night's  collection  session  lasted 
about  six  hours,  from  11:00  p.m  to  5:00  a.m.  local  time.  This  time  window  was  selected 
because  it  was  the  period  of  maximum  satellite  visibility  for  the  week  of  the  experiment. 
This  would  permit  the  tracking  of  at  least  four  satellites  simultaneously,  a  crucial  con- 
sideration from  the  standpoint  of  kinematic  operations  (at  least  four  satellites  are  needed 
to  solve  for  a  three-  dimensional  position  and  clock  bias).  Fortunately,  this  was  also  the 
time  when  ionospheric  effects  on  the  measurements  were  minimal.  All  GPS  measure- 
ments were  recorded  at  a  rate  of  one  per  second. 

The  MBPPE  static  GPS  positioning  data  was  collected  at  the  Naval  Postgraduate 
School's  Beach  Lab  in  Monterey,  California.  An  array  of  five  marks,  refered  to  as  the 
Doppler  Array,  was  established  to  serve  as  the  receiver  reference  points.  In  order  to  in- 
itialize the  kinematic  operations,  a  reference  array  (the  Lobos  Array)  was  established  at 
the  Monterey  Bay  Aquarium  Research  Institute's  pier  in  Moss  Landing,  California.  One 
of  the  points  in  this  array  would  provide  additional  static  GPS  measurements  that  would 
be  utilized  in  this  thesis.  Another  site,  established  on  the  roof  of  the  NPS  Mapping, 
Charting,  and  Geodesy  building  (in  Monterey),  would  supply  static  GPS  measurements 
for  the  Selective  Availability  test.  See  Figure  9  for  a  map  of  the  experiment  area. 
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Figure  9.      Map  of  Experiment  Area 
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B.     EQUIPMENT 

1.     Acquisition  and  Familiarization 

The  four  models  of  GPS  receiver  used  in  the  experiment  were  the  TI  4100,  the 
Trimble  4000ST,  the  Ashtech  LD  XII,  and  the  Magnavox  MX  4200.  Two  rubidium 
oscillators  accompanied  the  TI  4100  receivers.  For  a  summary  of  some  of  the  features 
available  on  each,  see  the  following  tables  (the  Magnavox  is  excluded  since  it  was  not 
used  in  point  positioning). 

Table  1.     ASHTECH  RECEIVER  FEATURES 


Standard  Features 

Comments 

No.  of  Channels 

12  dual  frequency  chan- 
nels 

Internal  RAM  Capacity 

6  MB  -  internal 
datalogging 

Recording  Interval 

1  second  or  more 

Frequencies 

Dual  Frequency  L1/L2 

Measurements 

L1/L2  Carrier  Phase,  C/A 
Code  Pseudorange  on  LI 

Table  2.     TRIMBLE  RECEIVER  FEATURES 


Standard  Features 

Comments 

No.  of  Channels 

8  dual  frequency  channels 

Internal  RAM  Capacity 

1  MB  -  used  external 
datalogger  (PC) 

Recording  Interval 

1  second  or  more 

Frequencies 

Dual  Frequency  L1/L2 

Measurements 

L1/L2  Carrier  Phase,  C/A 
code  Pseudorange  on  LI 

Two  of  each  type  of  receiver  was  acquired  either  from  the  vendors  themselves 
or  from  some  other  source.  The  Ashtech,  Trimble,  and  Magnavox  models  were  bor- 
rowed from  the  respective  vendors,  while  the  TI  4100s  were  on  loan  from  the  Texas 
Department  of  Highways  and  Public  Transportation.  The  receivers  were  obtained  a  few 
weeks  prior  to  the  data  collection  target  dates  so  that  the  experiment  participants  could 
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Table  3.     TI  4100  RECEIVER  FEATURES 


Standard  Features 

Comments 

No.  of  Channels 

4  dual  frequency  channels 

Internal  RAM  Capacity 

external  datalogging  via 
PC  only 

Recording  Interval 

1  second  or  more 

Frequencies 

Dual  Frequency  L1/L2 

Measurements 

L1/L2  Carrier  Phase,  P 

code  Pseudoranges  on 

L1/L2 

become  properly  familiarized  with  the  equipment.  One  person  was  assigned  the  respon- 
sibility of  learning  the  operational  aspects  of  a  particular  receiver  (as  applicable  to  ex- 
periment requirements).  Once  acquainted  with  the  equipment,  a  document  describing  the 
essential  details  of  the  receiver  operation  was  drafted  by  each  individual.  Then  all  other 
participants  were  cross-trained  on  the  different  receivers.  This  would  ensure  that  every- 
one would  be  able  to  start-up  operations  and  trouble  shoot  if  the  situation  arose. 
2.     Equipment  Set-Up  and  Data  Collection 

One  of  each  type  of  receiver  was  located  at  a  static  shore  site,  the  Naval  Post- 
graduate School's  Beach  Lab.  The  others  were  located  aboard  the  Research  Vessel 
POINT  SUR  where  the  kinematic  segment  of  the  experiment  was  conducted. 
a.     Beach  Lab 

Before  the  commencement  of  data  collection,  the  locations  of  the  stations 
over  which  the  individual  receiver  antennas  would  be  set-up  had  to  be  established.  At 
the  Beach  Lab  shore  site,  the  absolute  position  of  a  pre-existing  mark  (identified  as 
DOP)  had  previously  been  established  by  the  Transit  Doppler  method.  This  was  used 
as  the  reference  mark  from  which  the  positions  of  five  new  marks  would  be  established 
(designated  DOP1  -  DOP5).  This  configuration  of  closely  spaced  shore  marks  (about  five 
meters  separation  between  marks)  was  refered  to  as  the  Doppler  Array.  The  monuments 
were  set  and  the  positions  of  the  marks  were  determined  by  employing  conventional 
terrestrial  survey  techniques  using  a  steel  measuring  tape  to  to  obtain  distances  and  a 
Wild  T2000  theodolite  to  observe  horizontal  and  vertical  angles. 
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After  the  conclusion  of  the  experiment  data  collection  stage,  the  Transit 
Doppler  station  position  for  the  mark  DOP3  was  resurveyed  using  a  MX  1502  Transit 
Doppler  receiver  furnished  by  the  DMA.  Because  the  position  for  the  mark  had  been 
established  prior  to  the  October  1989  Loma  Prieta  earthquake,  a  post-earthquake  posi- 
tion for  the  mark  was  required.  The  Transit  Doppler  data  was  collected  over  four  days 
and  was  submitted  to  DMA  for  the  determination  of  the  DOP3  position  solution.  The 
positions  of  DOP3  and  the  other  stations  comprising  the  shore  site  array  were  then  up- 
dated based  on  the  new  information  produced  by  this  Transit  Doppler  survey.  These 
were  the  independent  position  solutions  to  which  the  GPS  solutions  were  compared. 

Each  receiver-antenna  was  assigned  to  an  individual  monument.  The 
Ashtech  antenna  was  set  up  over  DOP1,  the  Trimble  over  DOP2,  the  Magnavox  over 
DOP4,  and  the  TI4100  over  DOP5.  Except  for  the  second  night  of  data  collection,  when 
two  antenna  cables  were  inadvertantly  attached  to  the  wrong  receivers  (affecting  the 
Ashtech  and  TI  receivers),  this  configuration  remained  intact  throughout  the  exper- 
iment. See  Figure  10  for  the  Doppler  Array  layout. 
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Figure  10.      The  Doppler  Array:    The  NPS  Beach  Lab  Receiver  Monument  Sites, 
Ashtech  Set-up  Over  DOP1,  Trimble  Over  DOP2,  TI  4100  Over  DOP5 
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The  Krupp  Atlas  Polartrack  range  and  azimuth  laser  positioning  system  was  set  up  over 
mark  DOP3.  This  was  used  to  determine  the  reference  trajectory  of  the  ship  in  the 
kinematic  operations. 

The  heights  of  the  receiver  antennas  from  their  marks  were  measured  both 
at  the  start  and  end  of  nightly  operations.  These  measurements  are  necessary  to  adjust 
the  solution  from  the  antenna  electrical  center  to  the  mark. 

Meteorological  data  (i.e.,  temperature,  pressure,  and  relative  humidity)  was 
not  recorded  on  site  but  was  obtained  from  the  NPS  Meteorology  Department  for  the 
nights  of  the  experiment.  Because  the  location  of  the  School's  meteorology  recording 
station  is  in  close  proximity  to  the  Beach  Lab  site  (less  than  1000  meters),  the  weather 
data  should  closely  reflect  conditions  at  the  Beach  Lab. 
b.     The  Lobos  Site 

In  order  to  initialize  the  kinematic  operations,  the  pre-departure  and  post- 
arrival  positions  of  the  POINT  SUR,  GPS  antenna  array  (located  atop  the  mast  in  the 
crows  nest)  had  to  be  established.  The  positions  of  four  reference  points  located  in  the 
vicinity  of  the  POINT  SUR's  dock  were  determined  for  this  purpose.  This  station  group 
is  refered  to  as  the  Lobos  array  (See  Figure  11). 
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Figure  11.      The  Lobos  Array:     The  Moss  Landing  Monument  Sites,  TI  4100  Re- 
ceiver Set-up  Over  LOBOS3 

Following  the  experiment,  the  position  of  one  of  the  points,  LOBOS3,  was 
resurveyed  with  both  GPS  and  Transit  Doppler.  The  TI  4100  receiver  was  used  for  the 
GPS  survey.  The  GPS  data  was  collected  continously  over  a  three  day  period  (February 
5-7,  1991)  at  a  measurement  collection  rate  of  30  seconds.  This  supplied  three  inde- 
pendent data  sets,  one  for  each  day  of  the  survey.  The  Transit  Doppler  survey,  con- 
ducted from  January  31  -  February  4,  1991,  established  the  Transit  Doppler  position  for 
the  mark.  The  independent  GPS  solutions  are  compared  to  the  Transit  Doppler  solution 
in  later  sections. 
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The  height  of  the  TI  4100  antenna  above  the  mark  was  measured  at  at  the 
beginning  and  at  the  completion  of  the  three  day  collection  session.  Weather  data  was 
again  obtained  from  the  NPS  Meterology  Station.  Although  the  degree  to  which  this 
data  accurately  represents  conditions  at  Moss  Landing  is  not  known,  it  was  felt  it  should 
adequately  reflect  LOBOS3  weather  conditions. 

In  determining  the  point  positions  for  the  Beach  Lab  and  Moss  Landing 
sites,  the  precise  satellite  ephemerides  and  clock  models  were  obtained  from  DMA  for 
the  time  of  applicability.  This  provided  a  comparison  of  the  solutions  estimated  with  this 
information  versus  the  solutions  estimated  with  the  broadcast  ephemeris  and  clock  data 
collected  during  the  surveys. 
c.     Building  224 

As  of  early  July  1991,  Selective  Availability  (SA)  had  been  activated  by  the 
DoD.  GPS  data  was  collected  for  two  days  (July  16  and  17)  at  a  site  that  had  been 
previouly  established  on  the  roof  of  the  NPS  Mapping,  Charting,  and  Geodesy  center 
(site:  BLDG  224.3).  The  data  was  collected  with  two  different  TI  4100  receivers  at  a 
measurement  rate  of  30  seconds.  The  antenna  height  above  the  mark  was  recorded.  No 
weather  data  was  obtained  for  this  test.  The  data  was  then  processed  with  GASP  using 
the  broadcast  ephemeris  only  and  the  site  position  solution  was  compared  to  the  known 
site  position  in  order  to  evalute  the  effect  of  SA  on  the  computed  position. 
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III.     DATA  PROCESSING  WITH  GASP 

A.     PROCESSING  FLOW  OVERVIEW 

Before  the  raw  receiver  data  could  be  processed,  some  preliminary  operations  had 
to  be  performed  to  convert  the  data  into  a  format  suitable  to  the  GASP  program.  The 
raw  data  was  translated  on  a  PC  into  a  GASP  compatible  ASCII  format  (i.e.,  FICA  or 
RINEX).  The  actual  procedures  and  programs  involved  depend  upon  the  particular  re- 
ceiver and  are  covered  in  more  detail  in  the  next  section.  Once  the  data  had  been  con- 
verted and  all  other  required  data  assembled  in  the  proper  formats,  the  processing  was 
initiated. 

GASP  consists  of  two  major  program  units.  The  first,  called  the  GeoSTAR 
PREProcessor  or  STARPREP,  accepts  input  from  a  FICA  or  RINEX  data  file  and  from 
files  containing  satellite  ephemeris  and  clock  information,  meteorological  data,  and  sta- 
tion information  (such  as  antenna  height  above  mark  and  a  priori  station  coordinates). 
STARPREP  then  computes  and  applies  a  series  of  measurement  corrections  to  the  LI 
phase  biased  pseudorange  and  pseudorange  measurements.  Once  corrected,  these  meas- 
urements are  referred  to  as  two-frequency  corrected  measurements.  Also,  measurement 
time  tags  are  adjusted  to  reflect  the  correct  GPS  time  of  signal  transmission  and  the 
satellite  positions  are  interpolated  to  correspond  to  the  corrected  time  tags.  This  infor- 
mation, along  with  the  corrected  observations  and  the  values  computed  for  the  meas- 
urement corrections,  is  output  to  a  "point"  file.  A  station  information  file  is  also 
generated  containing  the  a  priori  station  position,  antenna  height,  and  other  relevant 
data.  These  files  are  subsequently  used  as  the  input  to  the  second  major  program  unit, 
GASP. 

The  GASP  unit  accepts  input  from  the  point  and  station  files,  forms  the  GASP  ob- 
servable, and  performs  the  estimation  for  the  point  position  in  CT  coordinates.  See 
Figure  12  for  a  general  depiction  of  the  processing  flow. 


36 


JZ]  ^y       FICA  o*"       i 

/T\  I     RINEX   Data      7 


C        Station  and       ^\ 
Mataoroloaical  J  - 


I 


STARPREP 


CSatanit*  EDhwTwktoa     \ 
and  OocK  Stataa         1 


c 


Point   File 


) 


i 


GASP 


I 


Point   Position 

GASP   Processing   Flow 


Figure  12.      GASP  Processing  Flow 

B.     PERSONAL  COMPUTER  (PC)  OPERATIONS 

The  operations  required  to  convert  the  receiver  data  into  a  GASP  acceptable  format 
varied  for  each  receiver.  An  outline  of  the  general  procedure  for  each  receiver  is  pre- 
sented. 

1.     TI  4100  Receiver  Data  Conversion 

The  raw  tracking  data  collected  with  the  TI4100  Basic  External  Processor  Pro- 
gram (or  BEPP)  operating  system  had  to  be  converted  to  the  Floating  point,  Integer, 
Character,  ASCII  (FICA)  format.  This  involved  the  use  of  utility  programs  developed 
by  the  University  of  Texas,  Applied  Research  Laboratory  (UT/ARL).  The  raw  tracking 
data  was  converted  to  the  Floating  point,  Integer,  Character  (FIC)  binary  format  by  the 
program  GS2FIC.  This  format  was  then  converted  to  the  FICA  format  by  the  FICFICA 
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program.  Because  the  Beach  Lab  data  was  collected  at  a  one  second  measurement  rate, 
it  was  necessary  to  decimate  the  data  to  30  seconds  in  order  to  alleviate  file  storage 
problems.  Also,  GASP  documentation  suggests  a  30  second  measurement  rate  and  over 
our  collection  periods  it  would  provide  more  than  a  sufficient  number  of  observations 
to  produce  meter  level  position  solutions. 

2.  Trimble  Receiver  Data  Conversion 

The  first  step  in  the  conversion  of  the  Trimble  raw  tracking  data  to  the  RINEX 
format  was  to  download  the  data  collected  by  the  receiver  into  a  restructured  binary 
format.  TRIMVEC  (Trimble  supplied  processing  software)  accomplished  this  task.  Then 
the  TRRINEX  programs,  developed  by  the  creators  of  the  RINEX  format  (a  team  from 
the  University  of  Bern,  Switzerland),  were  utilized  to  convert  the  Trimble  binary  format 
to  the  RINEX  format.  The  Beach  Lab,  Trimble  data  was  also  decimated  to  30  seconds. 
One  difference  between  the  FICA  and  RINEX  formats  should  be  noted  at  this  point,  the 
broadcast  ephemeris  and  clock  information  obtained  by  the  receivers  from  the  Naviga- 
tion Message  is  presented  as  separate  data  blocks  in  a  single  FICA  data  file.  In  the 
RINEX  format,  this  information  is  presented  as  a  separate  file  altogether.  This  fact  be- 
came important  when  considering  the  RINEX  modifications  to  GASP. 

3.  Ashtech  Receiver  Data  Conversion 

The  Ashtech  raw  data  also  had  to  be  downloaded  before  RINEX  conversion 
could  be  implemented.  The  Ashtech  GPPS  processing  software  performed  both  the 
downloading  and  the  RINEX  conversion.  ASHTORIN  was  the  program  used  to  per- 
form the  conversion.  It  too  permitted  the  Ashtech  Beach  Lab  measurements  to  be 
decimated  to  30  seconds. 

4.  PC  to  VAX  Data  Transfer 

The  observation  and  broadcast  ephemeris  data,  now  in  either  the  FICA  or 
RINEX  format,  had  to  be  transferred  to  the  NPS  Digital  VAX  computer.  This  is  where 
the  GASP  programs  and  program  code  resided  and  where  processing  would  be  done. 
The  data  was  transferred  to  the  VAX  via  9-track  tape  using  the  Overland  Data  Tape 
Software,  DEPOT  program.  The  precise  ephemerides  and  clock  information,  furnished 
by  the  DMA  for  the  weeks  of  applicability,  were  also  transfered  to  the  VAX. 
Meteorological  and  station  data  files  were  created  in  the  GASP  specified  formats  for  the 
different  days  and  receiver  stations.  See  Figures  13  and  14  for  examples  of  the 
meteorological  and  station  files. 
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Field  1  -  2-digit  year  identification 

Field  2  -  GPS  week  number 

Field  3  -  GPS  time  of  week 

Field  A  -  Barometric  pressure  in  millibars 

Field  5  -  Temperature  in  Celsius 

Field  6  -  Relative  humidity  in  percent 


Figure  13.       Sample  STARPREP  Meteorological  File  Contents  and  Record  De- 
scription 
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Record 

. 

Field  1  - 

Site  l.D.  (format  18) 

Field  2  - 

Record  type  (format  AI) 

freld  3  - 

Station  l.D.  (format  F10.0) 

Field  A  - 

Station  type  (format  Al) 

Field  5  - 

Station  name  (format  A30) 

Field  6  - 

Station  Latitude  in  rads  (format  DIG.  10) 

Field  7  - 

Station  Longitude  in  rads  (format  DIG.  10) 

Field  8  - 

Station  Latitude  in  degrees  (format  D15.9) 

Field  9  - 

Station  Longitude  in  degrees  (format  D15.9) 

Field  10 

-  Station  height  in  meters  (format  F10.6) 

Record  1 

•■  • 

Field  1  - 

Site  l.D.  2  (format  18) 

Field  2  - 

Record  type  (format  A  I) 

Field  3  - 

Site  type  (format  Al) 

Field  4  - 

Antenna  offset  north  of  mark  (format  F8.6) 

Field  5  - 

Antenna  offset  cast  of  mark  (format  F8.6) 

Field  6  - 

Antenna  height  offset  (format  F8.6) 

Field  7  - 

Year,  day  of  antenna  set  up  (format  15) 

Field  8  - 

Elevation  angle  cutoff  for  16  sv's  (format  16F4.1) 

Figure  14.  Sample  STARPREP  Station  File  Contents  and  Record 
Description:  Sv's,  referred  to  in  record  2,  field  8,  is  another  term  for 
satellites. 

C.     DATA  FILE  EDITS  AND  RUNSTREAM  GENERATION 

Prior  to  STARPREP  processing,  the  following  information  contained  in  the  FICA 
or  RINEX  observation  files  was  validated: 

•  the  a  priori  station  coordinates 

•  the  antenna  height 

•  the  year,  month,  and  day  of  session 

•  week-crossovers  , 
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•  the  PRN  numbers  of  the  satellites  tracked 

•  the  measurement  record  interval  (or  rate) 

Special  attention  had  to  be  paid  to  the  last  item  since  the  decimation  process  did  not 
automatically  update  the  measurement  rate  in  the  header  block  for  the  TI  4100  data. 

Much  of  the  data  collected  at  the  LOBOS3  site,  corresponding  to  periods  when 
fewer  than  three  satellites  were  available,  was  discarded.  This  still  provided  three  inde- 
pendent data  sets  (one  for  each  day  of  site  occupation)  of  between  eight  and  1 1  hours 
duration  at  a  30  second  recording  rate.  The  BLDG  224.3  data  set,  collected  contin- 
uously over  two  days  at  30  seconds,  was  divided  into  five  independent  subsets  ranging 
from  about  four  to  nine  hours.  The  data  sets  and  subsets,  time  spans,  and  the  number 
of  satellites  used  from  the  individual  data  sets  are  identified  in  the  following  tables. 


Table  4. 

BEACH    LAB    DATA    SET    IDENTIFICATION    AND 
TIME  SPANS 

Day 

Receiver 

Time  Span 
(hours) 

No.  of  Satellites 
Used 

338 

TI  4100 

6.2 

8 

Ashtech 

5.6 

9 

Trimble 

6 

9 

339 

TI  4100 

6.2 

8 

Ashtech 

6 

9 

Trimble 

6 

9 

340 

TI  4100 

7.1 

8 

Ashtech 

6.8 

9 

Trimble 

7.5 

9 

341 

TI  4100 

6.5 

8 

Ashtech 

6.9 

9 

Trimble 

6.7 

9 
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Table  5.     LOBOS3  DATA  SET  IDENTIFICATION 
AND  TIME  SPANS 


Day 

Time  Span  (hours) 

No.  of  Satellites 
Used 

36 

8.6 

8 

37 

10.8 

8 

38 

10.8 

8 

Table  6.     BLDG  224.3  DATA  SET  IDENTIFICATION  AND  TIME 
SPANS 


Day 

Subset 

Time  Span 
(hours) 

No.  of  Satellites 
Used 

197 

A 

4.8 

6 

B 

8.4 

6 

198 

C 

3.8 

5 

D 

7.0 

6 

E 

9.2 

6 

A  runstream  on  the  VAX  was  generated  that  signified  what  data  files  were  to  be 
used  as  input,  output  or  temporarily  created,  the  unit  numbers  assigned  to  each  file, 
which  of  the  available  satellites  were  to  be  used  as  the  source  of  observations,  whether 
broadcast  or  precise  ephemeris  and  clock  states  were  to  be  used,  the  corrections  to  be 
applied,  etc.  The  runstream  is  a  VAX  batch  file  that  begins  execution  of  the  STARPREP 
routine,  deletes  temporary  files  at  the  completion  of  processing  and,  in  general,  controls 
the  processing  operations.   See  Figure  15  for  a  sample  precise  ephemeris  runstream. 
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Figure  15.      Sample  STARPREP  Runstream  for  Precise  Ephemeris 
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D.    STARPREP 

1.  Input 

Much  of  the  information  presented  in  the  following  sections  on  STARPREP 
and  GASP  processing  is  found  in  Malys,  et  al  [Ref.  1]. 

STARPREP  accepts  input  data  from  an  observation  file  (pseudorange  and  car- 
rier phase  measurements  in  either  FICA  or  RINEX  format),  broadcast  or  precise 
ephemeris  and  clock  data  from  file  input  (or  in  the  case  of  the  FICA  broadcast 
ephemeris  data  as  part  of  the  FICA  observation  file),  a  meteorological  file,  and  a  station 
file.  The  carrier  beat  phase  measurements  are  converted  to  kilometers  by  multiplying 
by  the  nominal  LI  or  L2  wavelength.  In  the  TI  4100  FICA  observation  file,  the  time  of 
signal  transit  between  satellite  and  receiver  is  given  rather  than  the  actual  pseudorange. 
This  is  converted  to  the  pseudorange  by  multiplying  by  the  speed  of  light  in  a  vacuum. 
RINEX  presents  this  as  the  pseudorange  originally  so  that  no  conversion  to  distance 
units  is  necessary. 

2.  Error  Models 

The  time  and  data  corrections  applied  during  the  course  of  STARPREP  proc- 
essing will  be  briefly  discussed  and  the  error  models  used  to  compute  the  corrections 
presented.  This  material  is  originally  presented  by  Malys,  et  al  [Ref.  2:  pages  489-491]. 

a.  Time  of  Transmission 

Measurement  time  tags  are  adjusted  from  time  of  reception  (tr)  to  time  of 
transmission  (tx)  by, 

where,  p,  is  the  pseudorange  observation  at  /,  and 

c,  the  speed  of  light  in  a  vacuum,  equals  299792.458  kilometers  per  second.    See 
Appendix  C  for  a  list  of  the  constants  used  in  the  GASP  programs. 

b.  Satellite  Clock 

For  the  broadcast  ephemeris,  the  predicted  satellite  clock  parameters  avail- 
able in  the  navigation  message  are  used  to  compute  the  satellite  clock  offset  from  GPS 
time  by, 

tx  =  Oq  +  ax  (/,-/<>)  +  a2  {tx-t0)2 

where,  t0  is  the  time  of  applicability  of  %  au  a^ 
Op  is  the  predicted  satellite  clock  time  offset 
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ax  is  the  predicted  frequency  offset 
a2  is  the  predicted  frequency  drift. 

For  the  precise  ephemeris,  the  precise  clock  states  are  used  to  compute  the  satellite  clock 
correction  at  the  observation  epochs. 

c.  Receiver  Frequency  Offset 

The  TI  4100  receiver  LI  and  L2  reference  signals  are  offset  from  the  LI  and 
L2  carrier  frequencies  by  -6000  Hz  and  +  7600  Hz  respectively.  A  correction  is  computed 
which  removes  the  number  of  cycles  in  the  data  due  to  these  biases.  The  epoch  of  the 
initial  recorded  carrier  beat  phase  measurement  for  a  satellite  is  used  as  the  reference  for 
subsequent  carrier  phase  measurements.  For  example,  the  correction  for  LI  phase  data 
in  cycles  is  given  by, 

A^  — 6000(f-f,) 

where,  Nne  is  the  number  of  cycle  counts 

/  is  some  affected  epoch,  in  seconds,  and 

fj  is  the  initial  epoch  recorded  for  the  satellite,  in  seconds. 

The  Trimble  and  Ashtech  receivers  do  not  have  offsets  in  their  reference  signals  so  the 
values  assigned  to  the  offsets  were  set  equal  to  zero  for  RINEX  processing  runs  with 
these  receivers. 

d.  Ionosphere 

Signal  interaction  with  the  free  electrons  found  in  the  ionosphere  produces 
a  change  in  signal  path  length.  This  ionospheric  refraction  index  is  frequency  dependent 
so  a  comparison  of  measurements  on  LI  and  L2  may  be  used  to  derive  a  dual  frequency 
correction.  Defining  the  quantity  T  as 


r  = 


fui 


a  dual  frequency  ionospheric  correction  for  the  pseudorange  is  given  by, 

(Pru-PrL2) 


dp 


1i 
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This  can  only  be  applied  to  data  collected  with  the  TI  4100  receiver  since  it  is  the  only 
model  that  supplies  pseudoranges  from  the  L2  signal. 

The  dual  frequency  ionospheric  correction  for  the  carrier  phase  for  the  TI 
4100  receiver  is  given  by, 

dp*u  =  p*L1rp*u  +  y  Uu  ( -  6000)  -  xL1  (76oo):  m 

where,  Xx  and  X2  are  the  transmitted  wavelengths  and  At  is  the  interval  between  the  initial 
epoch  of  observation  and  the  epoch  being  corrected  [Ref.  2  :  p.  489].  For  the  Trimble 
and  Ashtech  receivers  with  no  receiver  frequency  offset,  this  reduces  to, 


*P*u  = 


yp<t>u    p<t>u) 


e.     Troposphere 

Tropospheric  effects  are  a  result  of  refraction  in  the  neutral  atmosphere  and 
are  not  frequency  dependent.  This  effect  may  be  separated  into  two  components:  the  dry 
component,  and  the  wet  component.  The  dry  component  comprises  about  90%  of  the 
total  effect  and  is  a  function  of  the  surface  atmospheric  pressure  and  satellite  elevation 
angle.  Estimation  of  the  wet  component  is  more  difficult  than  estimating  the  dry  com- 
ponent. It  depends  on  the  total  water  vapor  content  along  the  signal  path  and  hence 
on  the  temperature,  pressure,  and  humidity. 

In  GASP,  one  of  two  tropospheric  models  may  be  selected  to  compute  the 
integrated  tropospheric  correction,  the  Chao  or  Hopfield  models.  Both  models  use  sur- 
face weather  data  and  the  satellite's  elevation  angles  to  compute  the  correction.  The 
temperature,  pressure,  and  relative  humidity  are  used  to  compute  the  wet  and  dry  com- 
ponents of  the  zenith  tropospheric  delay  values  (Z„„,  Z^).  The  satellite  elevation  angles 
are  then  used  to  compute  the  wet  and  dry  multipliers  (F^,  Fwtl ).  The  total  correction, 
for  either  model,  may  be  generally  given  by, 

"Ptrop  =  \Adry  ''dry   >   ^wet  'wet) 

See  Chao  and  Hopfield  [Refs.  10,  11]  for  more  detailed  explanations  of  these  two  models. 
The  Hopfield  model  was  selected  as  the  tropospheric  model  for  all  the  re- 
sults produced  in  this  study.  A  few  of  the  similiarities  and  differences  between  the  two 
models  should  be  noted.  Both  models  produce  similiar  zenith  values  for  the  wet  and  dry 
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components.  The  major  difference  is  the  dependency  of  the  correction  on  the  satellite 
elevation  angles.  The  Chao  model,  originally  developed  for  use  in  arid  locations,  is  not 
as  accurate  as  the  Hopfield  model  in  estimating  the  wet  tropospheric  component  off  the 
zenith  [Ref.  12]. 

f.     General  Relativity 

If  using  the  broadcast  ephemeris,  the  general  relativity  correction  is  com- 
puted as  a  function  of  the  broadcast  orbital  elements  and  is  given  by, 

rl/2 


dprel  -  -2.0  [  -£*£—  ]  (e)  (a,/2)  [  sin(£a)] 


where,  GM,  the  product  of  the  universal  gravitational  constant  and  the  Earth's  mass, 
equals  3.986005  x  lO^m3/  sec2,  and 

e  is  the  eccentricity  of  the  satellite  orbit 

a  is  the  semi-major  axis  of  the  satellite  orbit 

Ea  is  the  eccentric  anomaly  at  the  observation  epoch. 

For  a  more  thorough  description  of  the  broadcast  orbital  elements  see  Appendix  A. 
If  the  precise  ephemeris  is  used,  the  correction  is  given  by, 

.  -2.0  X.  V 

dPrel  = c 

where,  X  is  the  satellite  position  vector  and 

Kis  the  satellite  velocity  vector  at  the  observation  epoch.    See  [Ref.  13]  for  a  full 
explanation  of  this  effect. 

g.     Earth  Rotation 

To  account  for  the  fact  that  the  Earth  is  rotating  while  the  signal  is  traveling 
from  satellite  to  receiver,  the  Earth  rotation  correction  is  given  by, 

dp„  =  "J"  t(X2  ~  Xrl)  *l  ~  C*l  "  *rl)  X21 

where,  co,  the  WGS  84  value  for  the  Earth's  rotation  rate,  equals  7.2921151467  x  10~s 
radians  per  second, 

X  is  the  interpolated  satellite  position  vector 

Xr  is  the  receiver's  a  priori  position  vector 
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1,  2  are  the  vector  components  along  the  X  and  Y  axes  of  the  WGS  84  reference 
frame. 

h.     Satellite  Antenna  Offset 

The  precise  ephemeris  contains  the  positions  of  a  satellites  center  of  mass. 
To  adjust  this  position  to  the  electrical  center  of  the  satellites  transmitting  antenna,  this 
satellite  antenna  offset  is  approximated  as, 

dPsa  =  R'S 

where,  R  is  the  range  vector  from  the  a  priori  receiver  position  to  the  satellite  and 

5  is  the  scaled  vector  from  the  satellite's  center  of  mass  to  the  satellite  sub-point. 
S  is  given  by, 

S  =  SA0  (ex) 

where,  SA0  is  the  L-band  satellite  antenna  offset  in  the  nadir  direction  (0.88m)  and 
ex  is  the  unit  vector  from  the  center  of  mass  of  the  Earth  to  the  satellite. 

The  broadcast  ephemeris  gives  the  position  for  the  phase  center  of  the  sat- 
ellite transmitting  antenna  [Ref.  14:  p.  83].  This  would  make  application  of  this  adjust- 
ment unnecessary  when  using  the  broadcast  ephemeris.  The  GASP  program  does  not 
distinguish  between  the  precise  and  broadcast  ephemeris  in  this  regard.  This  results  in  a 
range  error  that  is  common  to  all  satellites  when  the  broadcast  ephemeris  is  used.  Be- 
cause this  error  is  common  to  all  satellites,  the  GASP  algorithm  will  effectively  remove 
this  error. 

The  data  sets  were  originally  processed  with  this  correction  applied  for  both 
ephemerides.  Two  of  the  previously  processed  data  sets  were  reprocessed  using  the 
broadcast  ephemeris,  this  time  without  the  satellite  antenna  correction  applied.  The 
maximum  difference  for  any  component  over  the  two  data  sets  was  on  the  order  of  two 
centimeters.  It  can  be  reasonably  concluded  that  the  GASP  differencing  scheme  has 
eliminated  this  source  of  error. 
3.     Output 

The  STARPREP  preprocessor  generates  three  main  files;  the  point  file,  the  sta- 
tion file,  and  the  output  file.  GASP  accepts  the  point  and  station  files  as  input.  The 
output  file  contains  plots  of  all  the  applied  data  corrections  for  each  satellite  over  time, 
plots  of  the  uncorrected  range  and  carrier  phase  measurements  over  time,  and  a  sum- 
mary of  file  information.  This  file  is  helpful  in  identifying  problems  that  may  have  arisen 
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during  the  course  of  processing.  The  point  file  contains  the  two-frequency  corrected 
observations,  the  values  computed  for  the  individual  corrections,  the  corrected  meas- 
urement time  tags,  and  the  interpolated  satellite  positions  (in  CT  coordinates)  corre- 
sponding to  the  corrected  time  tag.  The  station  file  contains  the  a  priori  station 
coordinates  and  the  height  of  the  antenna  above  the  mark.  See  Figure  16  for  a  general 
summary  of  STARPREP  processing. 
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Figure  16.      STARPREP  Processing 

E.     GASP 

1.     GASP  Runstream  Generation 

Prior  to  GASP  processing,  a  GASP  runstream  was  created.  This  runstream 
specifies  the  station  and  point  files  to  be  used  as  input,  assigns  a  filename  to  the  GASP 
output  file,  and  if  desired,  permits  the  selection  of  non-standard  GASP  processing 
options.  Standard  options  preset  certain  processing  values  or  features  but  these  may  be 
overridden  by  the  operator.  See  Table  7  for  a  list  of  the  available  processing  options  and 
standard  values  and  features. 
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Table  7.     GASP   RUNSTREAM   PROCESSING  OPTIONS  AND 
STANDARD  VALUES  OR  FEATURES 


Processing  Options 

Standard 
Values  or 
Features 

Plot  to  Printer 

P 

Elevation  Angle  Cutoff  (Degrees) 

15 

RMS  Screening  Multiplier 

3.00 

Pseudorange  Editing  Tolerance  (meters) 

5.00 

Estimate  Fourth  Parameter 

N 

Minimum  of 'N'  Satellites  Per  Epoch  Pah- 

2 

Use  Offset  to  A  Priori  Station  Position 

N 

Standard  Deviation  on  GASP  Observables  (cm) 

20.0 

Standard  Deviation  on  A  Priori  Position  Components 
(km) 

0.05 

PRN  Number  as  'Base'  Sat,  Default  is  00  for  Sequenc- 
ing 

00 

Number  of  Batch  Least  Squares  Iteration 

3 

Number  of  Sequential  Estimation  Iterations 

1 

2.    Automatic  Data  Editing 

After  the  point  and  station  files  have  been  entered  into  GASP,  two  kinds  of  data 
editing  are  performed  to  ensure  consistency  between  the  pseudoranges  and  phase  biased 
pseudoranges  observed  from  each  satellite.  First,  the  corrected  pseudoranges  over  two 
successive  epochs  are  subtracted  and  the  corrected  phase  biased  pseudoranges  corre- 
sponding to  the  same  epochs  are  also  subtracted.  If  these  "delta"  ranges  differ  by  more 
than  some  user  specified  tolerance,  the  GASP  observables  (i.e.,  the  two  carrier  based 
phase  biased  pseudoranges)  are  rejected  from  the  data  set.  This  is  referred  to  as 
pseudorange  editing.  Note  that  this  is  the  only  use  of  the  corrected  pseudoranges.  They 
serve  only  to  monitor  the  phase  biased  pseudoranges  for  outliers  and  are  not  used  in  the 
estimation  process. 

For  the  processing  runs  performed  for  this  study,  the  standard  pseudorange  ed- 
iting tolerance  of  five  meters  was  initially  selected.  After  a  few  runs  with  the  Trimble 
and  Ashtech  receiver  data,  it  became  apparent  that  a  significant  percentage  of  the  data 
was  being  rejected  from  the  data  set  (on  the  order  of  35%  or  more).  In  order  to  avoid 
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this,  the  pseudorange  editing  tolerance  was  relaxed  to  ten  meters.  This  held  down  the 
amount  of  data  rejected  to  under  20  percent  in  most  cases.  There  are  two  possible  ex- 
planations for  why  such  a  large  amount  of  data  was  rejected.  The  first,  and  most  prob- 
able, an  ionospheric  correction  could  not  be  computed  for  the  Trimble  and  Ashtech 
pseudorange  data  since  no  pseudorange  measurements  were  supplied  on  the  L2  signal 
for  these  receivers.  This  would  produce  some  discrepancy  between  the  corrected 
pseudorange  and  the  phase  biased  pseudorange  where  an  ionospheric  correction  was 
applied.  A  second  possible  explanation  is  that  the  pseudoranges  are  observed  from  the 
C/A-code  modulations  rather  than  the  more  precise  P-code  modulations.  The  increased 
measurement  noise  on  the  C/A-code  observations  may  also  have  contributed  to  the 
higher  rejection  figures. 

The  second  edit  test  performed  on  the  data  uses  the  Root-Mean-Square  (RMS) 
of  the  residuals  of  the  previous  iteration  as  the  rejection  criterion.  The  RMS  is  initialized 
before  the  first  iteration.  Any  phase  biased  pseudorange  observables  with  a  residual 
greater  than  three  times  the  RMS  of  the  previous  iteration  are  rejected  from  the  data  set. 
These  editing  schemes  are  executed  prior  to  the  formation  of  the  GASP  observable. 
3.    The  GASP  Observable 

The  between-epoch  single  difference  equation  is  used  to  form  what  is  known  as 
the  GASP  observable.  Two  consecutive  carrier  based  phase  biased  ranges  from  the  same 
satellite  are  differenced.  Recall  that  this  between-epoch  single  difference  equation  .was 
given  by, 

&P<t>  =  Sp-c  {5dT  -  dxp) 

This  between-epoch  difference  is  then  differenced  with  the  corresponding  between-epoch 
difference  from  another  satellite.  The  resulting  equation  is  given  by 

ASpf  =  Spl  -  dpi 

=  (dp*  -  5pp)  -c(S4~  St*  +  8r*p  -  <5t£) 

where  a  and  /?  signifies  the  two  satellites.   It  may  be  expressed  more  conveniently  as, 

ASpj,  =  ASp  —  c  {A5tx  —  AStp) 

Since  the  satellite  clock  errors  and  the  errors  associated  with  various  other  sources  have 
been  modeled,  the  equation  becomes, 


51 


ASp  =  ASp^  +  c  {ASrx  —ASrp) 

The  value  ASp  is  the  GASP  observable.    See  Figure  17  for  a  conceptual  representation 
of  the  GASP  observable. 
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Figure  17.      The  GASP  Observable:     Differencing  Over  Satellites  Over  Two  Con- 
secutive Measurement  Epochs 
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This  type  of  differencing  scheme  reduces  the  effects  of  errors  in  satellite  clocks 
and  orbits,  and  removes  receiver  clock  errors  and  integer  ambiguities.  For  a  four  chan- 
nel receiver  with  all  channels  simultaneously  tracking  satellites,  three  GASP  observables 
can  be  formed  for  a  given  epoch  pair.  At  least  two  satellites  per  epoch  pair  are  required 
to  form  an  observable  for  that  epoch  pair.  For  each  epoch  pair,  one  satellite  is  used  as 
the  reference  from  which  the  others  are  differenced.  The  reference  satellite  is  selected 
sequentially,  so  that  for  every  new  epoch  pair  processed,  the  next  higher  satellite  PRN 
number  is  used  as  the  reference.  The  selection  sequence  cycles  back  to  the  lowest  PRN 
number  once  the  list  of  tracked  satellites  is  exhausted.  Individual  epochs  are  used  only 
once  to  form  an  observable. 

Reference  satellite  sequencing  reduces  correlation  among  the  observables.  Pre- 
cision, the  RMS  of  the  residuals,  and  the  variance-covariance  matrices  of  the  estimated 
parameters  are  all  improved  using  this  type  of  satellite  selection  rather  than  selecting  the 
satellite  with  the  most  stable  clock  as  the  reference.  [Ref.  1:  p.  24] 

For  example,  if  the  LI  phase  biased  pseudoranges  are  represented  by  p+  ,  four 
satellites  simultaneously  tracked  signified  by  superscripts  a,  ft,  y,  and  r\  (in  order  of  as- 
cending PRN  number),  and  four  consecutive  epochs  represented  by  subscripts  1,  2,  3, 
and  4,  a  between-epoch  difference  can  be  formed  for  each  satellite  given  by 

a  a  c    a 

P\-p\  =  bP\u 

The  three  GASP  observables  formed  by  satellite  sequencing  for  this  first  epoch  pair  are 
given  by, 

^pL-^pi^^pf 


bp\x-bp\=tep^ 


«5pJ12-^pJ12  =  A(5p^ 
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Then  using  /?  as  the  reference  satellite  for  the  next  epoch  pair  at  3  and  4,  the  GASP 
observables  are  given  by, 

The  next  epoch  pair  at  5  and  6  would  use  satellite  y  as  the  reference  and  epoch  pair  7 
and  8  would  use  satellite  rj.  Reference  selection  would  then  return  to  satellite  a  for  epoch 
pair  9  and  10. 

4.    Batch  Least  Squares  and  Sequential  Estimation 

After  the  GASP  observables  have  been  formed,  a  least  squares  technique  is  used 
to  estimate  the  station  position  components  in  the  CT  coordinate  system.  These  are  the 
only  parameters  estimated.  A  clock  correction  parameter  is  not  estimated.  It  is  assumed 
that  the  modeled  satellite  clock  states  have  sufficiently  aligned  the  satellite  clocks  to  GPS 
time  and  that  the  receiver  clock  error  has  been  removed  by  differencing.  The  interpo- 
lated satellite  positions  are  held  fixed  in  the  estimation  providing  the  reference  frame  in 
which  estimation  takes  place.  After  three  iterations  of  a  batch  least  squares,  the  esti- 
mated parameters,  the  variance-covariance  matrix  for  the  parameters,  and  the  RMS  of 
the  residuals  provide  the  input  to  the  second  estimation  step  (for  more  on  the  least 
squares  method  see  Appendix  B). 

This  is  a  sequential  estimation  algorithm  based  on  a  Kalman  filtering  routine. 
The  Kalman  state  is  the  receiver  position  vector.  It  does  not  include  any  clock  states. 
Fundamentally,  a  Kalman  filter  updates  measurements  from  one  observation  epoch  to 
the  next.  This  allows  the  parameter  estimates  and  covariances  to  be  updated  at  each 
measurement  epoch.  The  batch  method  provides  a  measure  of  data  noise  to  the  se- 
quential processor.  The  final  RMS  of  the  residuals  from  the  batch  processor  is  taken  as 
the  variance  of  one  GASP  observable  processed  through  the  sequential  processor.  Since 
the  sequential  processor  updates  the  station  coordinates  for  each  new  observable  proc- 
essed, the  estimated  station  coordinates  can  be  plotted  as  a  function  of  time.  Plots  of 
the  covariances  and  convergence  in  CT  coordinates  as  a  function  of  time  may  also  sup- 
plied. 

For  the  processing  runs  performed  in  this  study,  a  comparison  was  made  be- 
tween the  final  batch  least  squares  component  estimates  and  the  final  component  esti- 
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mates  that  had  passed  through  the  sequential  processor.    The  magnitude  of  the  mean 
component  differences,  averaged  for  all  processing  runs  using  the  precise  ephemeris,  was 
about  15  centimeters. 
5.     Output 

The  most  important  information  contained  in  the  output  file  generated  by 
GASP  is  the  estimated  station  position.  This  is  presented  as  the  X,  Y,  and  Z  coordinates 
of  the  CT  coordinate  system  and  in  geodetic  coordinates;  latitude,  longitude,  and  height 
(4>,  X,  h)  relative  to  the  WGS-84  ellipsoid.  Estimates  for  the  uncertainty  in  the  coordi- 
nates is  obtained  from  the  variance-covariance  matrix  for  the  estimated  parameters.  A 
correlation  matrix  is  computed  that  expresses  the  linear  independence  among  the  esti- 
mated parameters  from  the  parameter  variance-covariance  matrix.  A  posteriori  standard 
deviations  for  the  estimated  station  coordinates  are  obtained  by  taking  the  square  root 
of  the  diagonal  elements  (i.e.,  the  variances)  of  the  variance-covariance  matrix.  These 
provide  estimates  for  the  precision  of  the  estimated  point  position  components.  The 
GASP  processing  sequence  is  graphically  depicted  in  Figure  18. 
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Figure  18.      GASP  Processing 

F.     PROGRAM  MODIFICATIONS 

In  order  for  GASP  to  accept  the  data  collected  with  the  Trimble  and  Ashtech  re- 
ceivers, it  was  necessary  to  modify  the  program  to  allow  the  introduction  of  data  pre- 
sented in  the  RINEX  format.2  This  meant  altering  some  existing  program  units  and 
creating  two  new  program  modules  for  STARPREP.  The  first  reads  the  RINEX  obser- 
vation file.  The  STARPREP  subroutines  BLK006  and  BLK124  that  read  the  FICA 
measurement  and  station  data  blocks  were  used  as  models  in  building  the  RINEX  sub- 
routine. The  second  module,  called  RINNAV,  reads  the  RINEX  navigation  file  that 
contains  the  broadcast  ephemeris  information.  It  was  modeled  after  the  STARPREP 


2  RINEX  data  may  be  presented  in  one  of  two  possible  forms;  the  original  RINEX  format 
and  the  RINEX  2  (version  2)  format.  GASP  processing  runs  utilized  data  presented  in  the  original 
RINEX  format.  Program  code  designed  to  process  data  in  the  RINEX  2  format,  although  present, 
is  untested.  [Ref.  15)  , 
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subroutine  BLK009  that  reads  the  FICA  data  blocks  containing  the  broadcast  ephemeris 
data. 

To  differentiate  between  FICA  format  and  RINEX  format  input,  a  file  naming 
convention  was  implemented.  The  program  path  to  read  either  format  is  keyed  to  the 
first  four  characters  of  the  input  observation  filename.  For  the  RINEX  data  collected 
with  the  Ashtech  receiver,  the  first  four  characters  of  the  filename  must  be  RINA.  For 
the  RINEX  Trimble  data,  the  first  three  characters  must  be  RIN  and  the  fourth  char- 
acter anything  other  than  A.  For  FICA  data  collected  with  the  TI  4100,  the  first  three 
characters  may  be  anything  except  RIN.  The  reason  for  distinguishing  between  the 
RINEX  Ashtech  and  Trimble  data  is  due  to  a  problem  that  occurred  during  processing. 
A  wavelength  factor  of  one  was  specified  for  the  LI  signal  and  a  wavelength  factor  of 
two  specified  for  the  L2  signal  in  the  RINEX  observation  file  for  both  the  Ashtech  and 
Trimble  receivers.  Using  these  wavelength  factors  with  the  Ashtech  receiver  caused  all 
the  observables  to  be  rejected  by  the  pseudorange  edit.  This  was  because  of  an  incom- 
patibility between  the  internal  receiver  software  and  the  receiver  hardware.  Apparently, 
the  receiver  software  had  not  been  updated  to  reflect  changes  in  the  hardware.  To  cir- 
cumvent this  problem,  the  L2  wavelength  factor  for  the  Ashtech  receiver  is  reset  to  one 
based  on  the  RINEX  filename.  Once  the  receiver  software  is  updated,  the  program 
statement  that  resets  the  wavelength  factor  will  need  to  be  deleted. 

Also  keyed  to  the  filenames  are  values  assigned  to  the  receiver  frequency  biases.  The 
TI  4100  receiver  incorporates  receiver  frequency  biases  of -6000  Hz  on  the  LI  signal  and 
7600  Hz  on  the  L2  signal.  The  Ashtech  and  Trimble  receivers  do  not  incorporate  these 
biases. 

Because  the  Trimble  and  Ashtech  receivers  possess  more  tracking  channels  than  the 
TI  4100  and  could  thus  supply  simultaneous  observations  from  more  satellites,  we  felt 
we  should  take  advantage  of  this  feature  by  modifying  GASP  to  accomodate  data  col- 
lected on  these  additional  channels.  The  Trimble  receiver  records  observations  on  eight 
dual  frequency  channels  and  the  Ashtech  receiver  on  twelve.  Because  increasing  the  di- 
mensions of  the  program  arrays  might  cause  space  problems  and  adversely  affect  proc- 
essing, a  cautious  approach  was  taken  and  storage  expanded  to  accomodate  eight 
channels  per  epoch.  Basically,  this  was  accomplished  by  simply  redimensioning  the  ar- 
rays that  stored  the  input  measurements,  two-frequency  corrected  observations,  time 
tags,  and  interpolated  satellite  positions. 

The  program  revisions  to  include  more  channels  led  to  modifications  to  expand  the 
maximum  number  of  satellites  that  could  supply  observations  over  the  course  of  a  col- 
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lection  session.  The  previous  limit  of  eight  was  changed  to  twelve.  Thus,  the  number 
of  satellites  allowed  per  collection  session  should  not  be  a  limiting  factor  to  the  potential 
benefits  of  using  more  receiver  channels.  If  twelve  satellites  per  session  were  to  be  uti- 
lized, allowances  for  additional  files  associated  with  the  extra  satellites  (and  the  RINEX 
broadcast  ephemeris  file)  would  be  necessary.  The  maximum  number  of  permissible 
program  files  was  reset  from  60  to  62  and  program  statements  assigning  unit  numbers 
to  particular  files  or  setting  limits  on  the  unit  numbers  for  certain  file  groups  (such  as 
the  precise  ephemeris  and  temporary  files)  were  respecified.  Specifically,  the  unit  number 
assigned  to  the  RINEX  broadcast  ephemeris  file  is  11,  the  upper  and  lower  limits  on  the 
unit  numbers  for  the  first  set  of  temporary  satellite  files  changed  to  12  and  23,  limits  on 
the  precise  ephmeris  files  changed  to  24  and  35,  limits  on  the  output  data  files  changed 
to  36  and  39,  and  the  limits  on  the  last  group  of  temporary  files  changed  to  51  and  62. 
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IV.     RESULTS  AND  ANALYSIS 

In  order  to  provide  some  measure  for  the  accuracy  of  the  computed  positions,  the 
known  station  position  components  are  differenced  from  the  GPS  CT  estimated  position 
components  for  the  various  receivers  and  collection  sessions.  The  known  station  posi- 
tions are  the  Transit  Doppler  positions.  This  supplies  the  AX,  AY,  and  AZ  offsets  from 
the  "true"  positions.  Component  differences  AX,  AY,and  AZ  are  converted  to  AE,  AN, 
and  AU  (in  local,  cartesian  east,  north,  and  up  coordinates)  because  this  is  a  more  fa- 
miliar and  easily  comprehended  reference  frame.  The  following  discussion  and  analysis 
will  concentrate  on  the  tables  and  plots  that  present  the  results  in  the  local  east,  north, 
and  up  coordinate  system.  For  the  reader  interested  in  the  results  presented  in  CT  co- 
ordinates (X,  Y,  and  Z),  see  the  tables  in  Appendix  D. 

The  difference  results  are  displayed  in  the  form  of  tables  and  associated  target  plots. 
In  the  target  plots,  the  Transit  positions  are  represented  by  the  origin  and  the  symbols 
depict  the  values  computed  for  the  component  differences.  The  GASP  estimated  formal 
error  (one  sigma)  standard  deviations  are  presented  as  the  ±  terms  in  the  difference  ta- 
bles or  as  error  bars  in  the  difference  plots. 

The  known  station  positions  are  the  Transit  Doppler  derived  positions.  Mark  BLDG 
224.3  was  established  by  a  GPS  relative  positioning  survey  using  the  Transit  Doppler 
mark  DOP3  as  the  reference.  See  Appendix  E  for  a  summary  of  the  estimated  positions 
and  estimated  standard  deviations  on  the  positions  for  both  the  Transit  Doppler  and 
GPS  point  determinations. 

Repeatability,  also  defined  as  precision,  yields  a  measure  for  the  consistency  of  a 
given  set  of  position  results.  The  repeatability  is  determined  by  calculating  the  mean  and 
the  standard  deviations  (i.e.,  the  observed  errors)  on  a  set  of  position  component  differ- 
ences. Two  types  of  repeatability  were  determined.  One  averages  over  days  by  receiver 
accentuating  the  differences  between  the  receivers.  The  other  averages  over  receiver  by 
day  accentuating  the  difference  between  days.  In  the  repeatability  target  plots,  the  mean 
differences  are  represented  by  the  target  symbols  and  the  standard  deviations  are  re- 
presented by  the  error  bars. 

A.     BEACH  LAB  RESULTS 

Table  8  presents  the  accuracy  results  of  the  Ashtech,  Trimble,  and  TI  4100  receivers 
for  the  Beach  Lab  collection  sessions  using  the  broadcast  ephemeris  for  the  satellite  po- 
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sitions.  It  can  be  seen  from  this  table  and  the  corresponding  target  plots  (figures  19  and 
20)  that  the  GASP  estimated  components  give  fairly  good  agreement  with  the  Transit 
Doppler  estimates.  The  average  Root-Sum-Square  (RSS)  of  the  component  differences 
(i.e.,  the  magnitude  of  the  difference  vector)  is  at  about  the  two  meter  level.  The  origin 
is  within  the  ensemble  of  the  position  differences  and  formal  errors.  The  estimates 
produced  with  the  TI  4100  receiver  were  slightly  better  overall  than  the  Trimble  and 
Ashtech  estimates  but  all  were  comparable.  The  Trimble  results  are  the  most  widely 
dispersed  and  the  Ashtech  results  slightly  less  scattered.  The  TI  4100  results  exhibited 
the  least  dispersion  of  the  three  receivers. 

It  also  appears  that,  between  the  differences  computed  in  the  north  and  east  di- 
rections, the  greatest  dispersion  is  produced  in  the  east,  i.e.,  the  estimated  east  compo- 
nents are  less  precise  than  the  north.  For  all  three  components,  the  north  component 
shows  the  least  dispersion  i.e.,  the  best  agreement  is  between  the  GASP  estimated  north 
component  and  the  Transit  Doppler  north  component.  Comparison  of  the  AH  and  AU 
component  differences  shows  more  dispersion  in  the  vertical  than  in  the  horizontal. 

A  lack  of  overlap  between  the  GASP  estimated  formal  errors,  due  to  the  scatter  of 
the  component  differences,  is  also  conspicuous.  This  suggests  that  the  GASP  formal  er- 
rors underestimate  the  actual  observed  errors.  The  GASP  formal  errors  appear  to 
underestimate  error  in  the  east  component  to  a  greater  degree  than  the  north  compo- 
nent. The  dispersion  of  the  vertical  versus  the  horizontal  differences  suggests,  overall, 
the  GASP  estimated  errors  were  better  at  representing  the  horizontal  errors  than  the 
vertical.  This  visual  interpretation  is  supported  by  the  values  computed  for  the  bottom 
line  of  the  table  which  shows  the  mean  component  differences  and  standard  deviations 
on  the  means  over  all  receivers  and  days.  The  numbers  indicate  that  the  east  and  up 
components  were  underestimated  by  a  factor  of  about  three  and  the  north  by  a  factor 
of  approximately  two. 

Also,  as  compared  to  the  differences  computed  for  the  other  three  days,  the  day  339 
results  seem  "off.  This  may  be  due  to  poorer  broadcast  ephemeris  satellite  positions  for 
this  day. 
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Table  8.  BEACH  LAB  COMPONENT  DIFFERENCES  (AE,  AN,  AU)  USING  BROADCAST 
EPHEMERIS:  Solution  Differences  Between  GPS  and  Transit  Doppler  Methods,  Sec- 
ond Term  is  GASP  Estimated  Formal  Error,  RSS  is  Magnitude  of  the  Difference  Vector 


Day  of 
Year  1990 

Receiver 

AE  (m) 

AN  (m) 

AU  (m) 

AH  (in) 

RSS  (in) 

338 

TI  4100 

-0.6  ±.6 

-1.4  ±.6 

0.3  ±.5 

1.5  ±.8 

1.6 

Ashtech 

0.3  ±.6 

0.1  ±.4 

-1.8  ±.5 

0.4  ±.7 

1.9 

Trimble 

-0.6  ±.6 

1.1  +.4 

-2.5  ±.5 

1.3  ±.7 

2.8 

339 

TI  4100 

-0.4  ±  .6 

-1.5  ±.6 

1.7  ±.5 

1.6  ±.8 

2.3 

Ashtech 

3.3  ±.7 

1.0+.5 

-3.5  ±  .5 

3.4 +  .8 

4.9 

Trimble 

2.9  ±.6 

1.2 +  .5 

-3. 9  ±.5 

3.1  ±.7 

5.0 

340 

TI  4100 

-0.6  ±  .5 

-0.6  ±.5 

0.5  ±  .5 

0.9  ±  .7 

1.0 

Ashtech 

0.8  ±.6 

-0.1  ±.5 

-1.9  ±.5 

0.8  ±.6 

2.1 

Trimble 

0.4  ±  .5 

0.8  ±  .4 

-1.5  ±.4 

0.8  ±  .6 

1.7 

341 

TI  4100 

1.4 +  .5 

-0.8  ±  .5 

-0.5  ±  .4 

1.6  ±.7 

1.7 

Ashtech 

-1.8  ±.5 

0.3  ±.4 

-0.5  ±  .5 

1.8 +  .6 

1.9 

Trimble 

0.3  ±  .5 

-0.0  ±  .4 

0.2 +  .5 

0.3  ±  .6 

0.4 

Mean  Component  Differ- 
ence and  Standard  Devi- 
ation about  Mean 

0.5  ±1.5 

0.0  ±  .9 

-1.1  ±1.7 

1.5  ±  1.0 

2.2  ±  1.3 
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Table  9  and  Figures  21  and  22  present  the  accuracy  results  of  the  three  receivers 
using  the  precise  ephemeris.  It  can  be  seen  that  the  position  estimates  produced  with  the 
precise  ephemeris  give  better  agreement  to  the  Transit  Doppler  positions  than  the 
broadcast  ephemeris.  Overall,  the  precise  ephemeris  results  are  more  accurate  and  pre- 
cise in  every  component  difference  as  demonstrated  by  comparison  of  the  mean  com- 
ponent differences  and  the  standard  deviations  computed  about  the  mean  differences. 
The  overall  RSS  is  at  about  the  1.5  meter  level.  Again,  the  origin  is  within  the  ensemble 
of  all  position  differences  and  formal  errors. 

The  estimates  produced  by  the  three  receivers  are  all  comparable,  with  the  Ashtech 
and  Trimble  giving  somewhat  better  results  than  the  TI  4100  receiver  over  the  four  days 
of  observations  especially  in  terms  of  difference  dispersion.  The  component  differences 
for  the  Ashtech  receiver  are  all  fairly  tightly  clustered  about  the  origin  while  the  Trimble 
and  TI  4100  receivers  show  greater  dispersion  especially  in  the  east  component.  The 
TI  4100  receiver  also  shows  the  greatest  dispersion  in  the  north  and  up  directions. 

As  seen  in  the  broadcast  ephemeris  results,  comparisons  of  dispersion  between  the 
north  and  the  east  reveal  that  dispersion  is  greater  in  the  east  than  in  the  north  direction. 
Also  conspicuous  is  the  greater  dispersion  in  the  vertical  than  in  the  horizontal. 

Again,  we  see  that  the  GASP  estimated  formal  errors  have  underestimated  the  ac- 
tual error  as  determined  by  the  dispersion  of  the  differences.  For  the  precise  ephemeris 
estimates,  GASP  formal  error  underestimates  the  actual  error  more  in  the  east  and  up 
components.  It  represents  the  error  in  the  north  component  fairly  well.  The  actual  error 
on  the  east  and  up  components  is  underestimated  by  a  factor  of  about  two.  The  re- 
lationship for  the  error  in  the  north  component  is  almost  one-to-one. 
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Table  9.     BEACH    LAB   COMPONENT   DIFFERENCES   (AE,   AN,   AU)    USING    PRECISE 

EPHEMERIS:     Solution  Differences  Between  GPS  and  Transit  Doppler  Methods,  Sec- 
ond Term  is  GASP  Estimated  Formal  Error,  RSS  is  Magnitude  of  the  Difference  Vector 


Day  of 
Year  1990 

Receiver 

AE  (m) 

AN  (in) 

AU  (m) 

AH(m) 

RSS(m) 

338 

TI  4100 

1.3  ±.6 

1.3  ±.6 

-1.6  +  .5 

1.9  ±.8 

2.4 

Ashtech 

-1.2 +  .6 

-0.6  ±  .4 

-0.7 +  .5 

1.4  ±.7 

1.5 

Trimble 

—  1.9  ±  .6 

0.3  ±  .4 

-0.8  ±  .4 

1.9  ±.6 

2.1 

339 

TI  4100 

-0.5  ±.6 

-0.1  ±.6 

-0.6  ±  .5 

0.5  ±.7 

0.7 

Ashtech 

-0.3  ±  .6 

-0.6  ±  .5 

-0.0  ±  .5 

0.7 +  .7 

0.7 

Trimble 

-0.0  ±.6 

-0.1  ±.4 

-1.4  ±.4 

0.1  +  .5 

1.4 

340 

TI  4100 

-I.6±.5 

-0.0  ±  .5 

— 1.6±  .5 

1.6  ±.6 

2.3 

Ashtech 

-0.0  ±  .6 

-0.4  ±  .5 

-2.3  ±.5 

0.4  ±  .5 

2.3 

Trimble 

-0.4  ±  .5 

0.2  ±  .4 

-1.3  ±.4 

0.4  ±  .6 

1.4 

341 

TI  4100 

2.1 +.6 

0.8  ±.5 

1.6  ±.5 

2.2  ±.7 

2.8 

Ashtech 

-0.2  ±  .6 

0.6  ±  .5 

-0.1  ±.5 

0.7  ±  .6 

0.7 

Trimble 

1.4  ±.6 

0.2  ±.4 

0.5  ±  .5 

1.4 +  .6 

1.5 

Mean  Component  Differ- 
ence and  Standard  Devi- 
ation about  Mean 

-0.1  ±1.2 

0.1  ±.6 

-0.7  ±1.1 

1.1  ±.7 

1.6  ±.7 

65 


< 

2 

0 

-2 

-4 

i     i      1      i      i      i      i      1      i      i      i      i      1      i      i      i      i      1      i      i" 

1      i      i      i      i 

l         c 

'      1 
i 

i     i      i     i 
i 

• 

1         - 

3               1- 
-    1 

« 

1 

1             - 

j       i 

i     i     i     i 

>— 

P-I 

1 

"        1 

|— 

i — <>= 

1 — CS 1            L- 

<M         I < 

"1*      1 

* 1 

,           1 

- J 

-1 — Sfc 

L      . 

-1      ' 

i 

1— 

iL      1 

1  1 — r±J* 

1 

1      .      .      ,      .      1      .      ,      . 

RECEIVER 
s     TI  4100 
x     Ashtech 
o     Trimble 

-2                        0 

AE  (m) 

2 

4 

Figure  21.      AEast  versus  ANorth  for  Beach  Lab,  Precise  Ephemeris  Data 


66 


<] 


i     i     i     i      1 

-i—  "i      i      i      1      i      i      i      i      |      i      i 

- 

- 

2 

1            e:            1 

-                            I            :  j            1 

i 

- 

if 

- 

0 

■!■-      n 

1 

i 

I 

3                1  . 

1                 P 
1 

*-+- 

— 1                                                          ' 
...   1 

it — it 

L 

2 

i — 

i 

j— 

■H 

i 

-  -L 

- 

- 

- 

- 

- 

4 

— 

— 

- 

-    ,     ,     .     ,     1 

i      i      i      i      1      i      i      i      i      1      i      i 

2  4 

AH  (m) 


RECEIVER 
25     TI  4100 
x     Ashtech 
o     Trimble 


Figure  22.      AHor  versus  A  Up  for  Beach  Lab,  Precise  Ephemeris  Data 


67 


Table  10  and  Figures  23  and  24  display  the  Beach  Lab  repeatability  results  that  av- 
erage over  days  by  receiver  for  both  the  precise  and  broadcast  ephemerides.  The  means 
and  standard  deviations  on  the  means  show  that  geodetic-quality  point  positions  are 
achieved  using  the  precise  ephemeris.  This  is  especially  true  for  the  results  produced  with 
the  Ashtech  and  Trimble  receivers.  The  standard  deviations  i.e.,  observed  errors  are 
typically  within  or  around  one  meter  on  all  components  for  these  receivers.  The  accura- 
cies, as  determined  by  the  RSS  of  the  mean  component  differences,  are  within  one  meter 
for  all  receivers  demonstrating  that  agreement  to  the  Transit  Doppler  estimates  are 
within  the  computed  noise  levels. 

The  TI  4100  precise  ephemeris  results  are  less  precise  than  either  the  Ashtech  or 
Trimble  precise  ephemeris  results  and,  interestingly,  are  less  precise  than  the  TI  4100 
broadcast  ephemeris  results.  The  reason  for  this  is  not  known.  However,  examination 
of  the  RSS  of  the  mean  component  differences  reveals  that  the  overall  accuracy  of  the 
TI  4100  precise  ephemeris  estimates  is  better  than  the  broadcast  ephemeris  estimates. 
This  is  presumably  due  to  increased  systematic  error  associated  with  the  broadcast 
ephemeris  satellite  positions  and  clock  states. 

The  broadcast  ephemeris  estimates  show  that  geodetic-quality  point  positions  are 
not  obtained  using  the  broadcast  satellite  positions  and  clock  states.  Higher  values 
computed  for  the  RSS  of  the  mean  differences  show  reduced  accuracy  compared  to  the 
precise  ephemeris  estimates.  Also,  the  broadcast  ephemeris  estimates  are  generally  more 
imprecise  than  the  precise  ephemeris  estimates.  As  previously  mentioned,  this  was  not 
true  for  the  TI  4100  receiver  however. 

It  is  also  clear  that  averaging  the  component  differences  over  days  has  significantly 
improved  the  accuracy  and  precision  for  all  receivers  for  both  ephemerides.  The  indi- 
vidual components  show  improvement  and,  as  a  consequence,  the  RSS  of  the  averaged 
components  is  also  improved.  This  demonstrates  the  importance  of  averaging  over  many 
independent  position  estimates  to  suppress  the  effects  of  random  error. 
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The  error  bars,  representing  the  actual  observed  error  (the  standard  deviations  about 
the  mean  component  differences),  show  significant  asymmetry.  This  dramatically  illus- 
trates the  reduced  precision  in  the  east  component  as  compared  to  the  north  component 
and  in  the  vertical  component  when  compared  to  the  horizontal.  It  also  illustrates  that 
the  GASP  estimated  formal  errors  substantially  underestimated  the  observed  or  true  er- 
rors. For  the  precise  ephemeris  estimates,  errors  on  the  north  component  for  all  three 
receivers  displayed  close  agreement  to  the  GASP  estimated  errors.  The  Ashtech  receiver 
showed  the  best  agreement  in  the  east  component  and  the  Trimble  receiver  in  the  up. 
The  broacast  ephemeris  estimates  also  displayed  close  agreement  in  the  north  for  all  re- 
ceivers, but  the  east  and  up  components  showed  even  poorer  agreement  than  the  precise 
ephemeris. 
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Table  10.     BEACH  LAB  REPEATABILITY  (AE,  AN,  AU):  AVERAGE  OVER  DAYS  BY 

CEIVER:     Mean  RSS  is  Magnitude  of  the  Mean  Difference  Vector 


Ephemeris 

Receiver 

Mean  AE 
(m) 

Mean  AN 
(m) 

Mean  AU 
(m) 

Mean  AH 
(m) 

Mean  I 
(m) 

Broadcast 

TI  4100 

0.0+1.0 

-1.1  ±.4 

0.5  +  .9 

1.1  ±.5 

1.2  i 

Ashtech 

0.7  +  2.1 

0.3 +  .5 

-1.9+1.2 

0.7  ±2.1 

2.1  ± 

Trimble 

0.8  ±1.5 

0.8  ±.6 

-1.9+1.8 

1.1  ±1.4 

2.2  ± 

Precise 

TI  4100 

0.3  +  1.7 

0.5  +  .7 

-0.5  ±1.5 

0.6  ±1.5 

0.8  ± 

Ashtech 

-0.5  ±  .5 

-0.2  ±  .6 

—  0.8  ±  1.1 

0.5  ±  .7 

0.9  ±] 

Trimble 

-0.2+1.3 

0.2  +  .2 

-0.7  ±  .9 

0.3  ±1.2 

0.8  ±| 
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Table  1 1  and  Figures  25  and  26  display  the  Beach  Lab  repeatability  results  that  av- 
erage over  receivers  by  day  for  both  the  precise  and  broadcast  ephemerides.  The  table 
and  the  plots  show  that  the  day  339  results,  produced  using  the  broadcast  ephemeris, 
were  the  most  inaccurate  and  imprecise  of  all  the  daily  estimates.  Overall,  the  broadcast 
ephemeris  results  showed  much  less  precision  than  the  precise  ephemeris  estimates.  The 
precision  on  the  positions  estimated  with  the  precise  ephemeris  was  roughly  two  times 
better  than  those  estimated  with  the  broadcast  ephemeris.  This  demonstrates  the  in- 
consistency of  the  broadcast  ephemeris  in  estimating  geodetic-quality  point  positions. 

Again,  the  plots  show  asymmetry  in  the  magnitudes  of  the  error  bars.  The  east 
component  displays  greater  imprecision  than  the  north  component.  Comparison  of  the 
vertical  and  horizontal  components  does  not  reveal  any  obvious  overall  asymmetry  ex- 
cept in  a  few  individual  cases.  This  also  shows  that,  again,  the  GASP  estimated  formal 
errors  underestimate  the  true  errors.  The  computed  values  demonstrate,  for  the  precise 
ephemeris,  the  north  and  up  formal  errors  best  estimate  the  observed  errors  for  all  days. 
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Table  11.     BEACH  LAB  REPEATABILITY  (AE,  AN,  AU):  AVERAGE  OVER  RECEIVERS! 

DAY:     Mean  RSS  is  Magnitude  of  the  Mean  Difference  Vector 

Ephemeris 

Day 

Mean  AE 
(m) 

Mean  AN 
(m) 

Mean  AU 
(m) 

Mean  AH 
(m) 

Mean  E 
(m) 

Broadcast 

338 

-0.3  ±.6 

-0.1  +  1.3 

-1.4+1.5 

0.3  ±.8 

1.4  ±1 

339 

1.9  ±2.0 

0.2+1.5 

-  1.9  +  3.1 

1.9  +  2.2 

2.7  +  2 

340 

0.2  ±.7 

-0.0  ±  .7 

-  1.0  ±  1.3 

0.2 +  .7 

1.0  ±1 

341 

-0.0+  1.6 

-0.1  ±.6 

-0.3  ±  .4 

0.2  ±.8 

0.3  ±. 

Precise 

338 

-0.6±  1.7 

0.4+1.0 

-1.0 +  .5 

0.7  ±  1.9 

1.3  ±1 

339 

-0.3  ±  .2 

-0.3  +  .3 

-0.7  +  .7 

0.4  +  .4 

0.8  ±. 

340 

-0.7  ±  .8 

-0.1  ±.3 

-1.8  ±.5 

0.7 +  .8 

1.9  ±. 

341 

1.1  +  1.2 

0.5  ±.3 

0.7  +  .9 

1.2+1.2 

1.4+1 
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Figure  26.      Mean  AHor  versus  Mean  AUp  for  Beach  Lab  Repeatability:  Average 
Over  Receivers  by  Day 
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B.     LOBOS3  RESULTS 

The  TI  4100  receiver  was  the  only  model  used  to  collect  data  at  the  LOBOS3  site. 
Independent  data  sets  were  created  for  each  of  the  three  days  of  site  occupation.  Results 
produced  by  the  broadcast  and  precise  ephemeris  are  compared  for  each  of  the  inde- 
pendent data  sets. 

Table  12  and  Figures  27  and  28  displaying  the  LOBOS3  component  differences  for 
both  ephemerides  clearly  demonstrate  the  advantage  of  using  the  precise  ephemeris  to 
determine  the  point  position.  Accuracy  is  improved  in  every  component.  The  broadcast 
ephemeris  estimates  show  much  lower  levels  of  accuracy  and  precision  and  there  appears 
to  be  a  bias  present  in  the  south-east  and  up  directions.  Errors  in  the  predicted  satellite 
positions  or  in  the  broadcast  satellite  clock  models  are  the  probable  cause  of  this  bias. 

The  precise  ephemeris  results  show  very  good  agreement  with  the  Transit  Doppler 
estimates  and  can  definitely  be  considered  geodetic-quality.  As  with  the  Beach  Lab  re- 
sults, the  precise  ephemeris  estimates  show  greater  dispersion  in  the  east  component 
than  in  the  north.  Any  difference  in  dispersion  between  the  vertical  and  the  horizontal 
does  not  appear  to  be  significant  for  either  ephemeris.  Also  noteworthy,  the  TI  4100 
gave  much  better  precise  ephemeris  results  at  this  site  than  at  the  Beach  Lab  and  the 
broadcast  ephemeris  results  are  much  worse  overall. 

Because  of  the  peculiar  TI  4100,  precise  ephemeris  results  obtained  at  the  Beach 
Lab,  it  would  be  desirable  to  compare  the  precise  ephemeris  results  obtained  with  the 
TI  4100  receiver  at  LOBOS3  to  those  obtained  by  the  other  two  receivers  at  the  Beach 
Lab.  Direct  comparison  of  the  results  requires  some  caution,  however,  because  of  the 
longer  lengths  of  the  data  sets  processed  for  LOBOS3.  With  this  caveat  in  mind,  com- 
parison of  the  LOBOS3  TI  4100  estimates  to  the  Ashtec  and  Trimble  Beach  Lab  esti- 
mates reveals  that  better  results,  in  terms  of  both  accuracy  and  precision,  were  produced 
with  the  TI  4100  receiver  for  all  components  at  LOBOS3.  This  may  be  due  to  better 
precise  ephemeris  satellite  positions  provided  for  this  particular  week.  Again,  it  should 
be  stressed,  although  the  LOBOS3  precise  ephemeris  results  with  the  TI  4100  were  best 
overall,  all  three  receivers  produced  geodetic-quality  point  positions  with  the  precise 
ephemeris. 
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Table  12.     LOBOS3  COMPONENT  DIFFERENCES  (AE,  AN,  AU):     Solution  Differences  I 
tween  GPS  and  Transit  Doppler  Methods  for  Broadcast  and  Precise  Ephemerides, 
ond  Term  is  GASP  Estimated  Formal  Error,  RSS  is  Magnitude  of  the  Difference  Vei 


Ephemeris 

Day  of 
Year  1991 

AE(m) 

AN  (m) 

AU(m) 

AH(m) 

RSS(i 

Broadcast 

36 

1.4  ±.6 

-1.6  +  .5 

1.7 +  .5 

2.1 +  .8 

2.7 

37 

3.4  ±.8 

-4.6  ±.6 

4.9  ±.6 

5.7  ±1.0 

7.5 

38 

1.3  ±.7 

-3.3 +  .5 

5.8 +  .5 

3.6 +  .7 

6.8 

Precise 

36 

-0.4  ±  .6 

-0.4  ±  .5 

0.6  ±.5 

0.6  ±.8 

0.9 

37 

0.4 +  .6 

-0.2  ±  .4 

-0.6  ±  .4 

0.5  ±  .8 

0.7 

38 

1.3  ±.6 

0.1  ±.5 

0.2  ±.5 

1.3  ±.7 

1.3 

78 


<3 


-4 


-6 


n — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — r 


I I 


f£& 


_J_ 


-2024 

AE  (m) 


PRECISE  EPHEMERIS 


DAY 

X 

36 

o 

37 

D 

38 

-2    — 


-4    — 


t — n — i — r- 1 — i — r— i — i— i — i — i    i    I — i    i    i    i    I    i 


-6 


I  j  I 


i 


\—h- 1 


-1—1 '    ■   '    i   ■   ■    ■    I I    ■    I    ■    ■    I 


-2024 

AE  (m) 


BROADCAST  EPHEMERIS 


DAY 

X 

36 

o 

37 

a 

36 

Figure  27.      AEast  versus  ANorth  for  LOBOS3 


vr 


79 


7T — i    i    i    I    i    i    i    i    I    I    I    i    i    I    i    i    I    I    I — r 


4    — 


-2 


-4 


-6  LJ-J 


_i i i i i_j i i   i   i   i   i i t 


2  4  6 

AH  (m) 


PRECISE  EPHEMERIS 


DAY 

X 

36 

o 

37 

D 

38 

I   I    I    I   I    I    I    I    I    I   I    I    I    I    I   I    I   I   I    I 

I— J— I 


-2 


-4    — 


* 


y—4 — i 


-6 


i 


I ♦ 1 


■   i   ■   ■   I   i   t   ■   i   I   i   ■   ■   ■   I   ■   ■■   ■   I   ■   ■ 


0  2  4  6  6 

AH  (m) 


BROADCAST  EPHEMERIS 


DAY 

X 

36 

o 

37 

o 

38 

Figure  28.      AHor  versus  AUp  for  LOBOS3 


80 


The  LOBOS3  repeatability  results  averaging  over  days,  displayed  in  Table  13  and 
Figures  29  and  30,  clearly  shows  the  inaccuracy  and  imprecision  of  the  broadcast 
ephemeris  results.  The  south-east-up  bias  alluded  to  earlier  can  also  be  clearly  seen  in 
the  broadcast  ephemeris  estimates.  Comparison  of  the  error  bars  on  the  averaged  precise 
ephemeris  estimates  shows  the  characteristic  east  versus  north  asymmetry.  No  signif- 
icant vertical  versus  horizontal  asymmetry  is  evident,  however.  The  magnitudes  of  the 
error  bars  on  the  precise  ephemeris  estimates  show  that,  in  this  case,  the  GASP  esti- 
mated formal  errors  more  closely  represented  the  true  errors.  However,  the  formal  error 
estimates  produced  for  the  broadcast  ephemeris  results  underestimated  the  true  errors 
by  a  factor  of  approximately  two  for  the  east  component,  three  for  the  north,  and  four 
in  the  vertical.  The  precise  ephemeris  statistics  demonstrate  that  by  averaging  over  in- 
dependent estimates,  the  accuracy  and  precision  can  be  significantly  improved. 

C.     REGIONAL  AGREEMENT  BETWEEN  GPS  AND  TRANSIT  DOPPLER 

In  order  to  establish  a  regional  reference  for  the  expected  agreement  between  GPS 
and  Transit  Doppler  position  estimates  in  the  Central  California  area,  the  overall  accu- 
racy and  precision  of  the  position  differences  using  the  precise  ephemeris  are  evaluated. 
The  precise  ephemeris  mean  component  differences  and  standard  deviations  over  all 
Beach  Lab  and  LOBOS3  data  sets  were  computed  to  provide  values  for  this  expected 
agreement.  These  values,  obtained  by  averaging  over  15  data  sets,  are  mean  AE 
=  0.21  +  1.1,  mean  AN  =  0.09  ±  .52,  and  mean  AU  =  -0.54  ±  1.03. 
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Table  13.     LOBOS3  REPEATABILITY  (AE,  AN,  AU)  AVERAGE  OVER  DAYS  FOR  TI  4 
RECEIVER:     Mean  RSS  is  Magnitude  of  the  Mean  Difference  Vector 


Ephemeris 

Mean  AE  (m) 

Mean  AN  (m) 

Mean  AU  (m) 

Mean  AH  (m) 

Mean  RS 
(m) 

Broadcast 

2.0  ±1.2 

-3.2  ±  1.5 

4.1  ±2.2 

3.8  ±1.9 

5.6  ±2.9 

Precise 

0.4  ±.8 

-0.2  ±.3 

0.1  ±.6 

0.5  ±  .8 

0.5  ±  .9 
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D.     POSITIONING  RESULTS  -  SELECTIVE  AVAILABILITY  ACTIVATED 

Table  14  and  Figures  31  and  32  present  the  solution  differences  between  the  coor- 
dinates previously  established  for  mark  BLDG  224.3  and  the  coordinates  estimated  for 
the  mark  using  the  broadcast  ephemeris  after  the  activation  of  Selective  Availability. 
From  comparison  of  the  broadcast  ephemeris  results  obtained  at  the  Beach  Lab  and 
LOBOS3  sites  to  the  results  displayed  here,  it  is  quite  obvious  that  the  activation  of  SA 
has  severely  impacted  the  position  estimate.  Both  the  accuracy  and  precision  of  the  re- 
sults have  been  adversely  affected. 

The  average  RSS  of  the  estimated  positions  over  the  five  data  sets  is  about  16  me- 
ters. When  contrasted  to  the  approximately  six  meter  average  RSS  computed  for  the 
LOBOS3  site  using  the  broadcast  ephemeris,  some  indication  of  the  magnitude  of  this 
effect  is  realized.  The  east  and  up  components  are  most  affected.  The  north  component 
seems  to  be  much  less  effected  by  SA  than  the  other  two  components.  Precise 
ephemerides  were  not  obtained  for  this  survey  period  so  a  comparison  of  the  broadcast 
ephemeris  results  to  the  precise  ephemeris  estimated  position  cannot  be  performed. 
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Table  14.  BLDG  224.3  COMPONENT  DIFFERENCES  (AE,  AN,  AU):  Selective  Availabi 
Test  -  GPS  Solution  Differences  Between  Previously  Determined  Coordinates  and  < 
ordinates  Determined  After  Activation  of  S/A  (Using  Broadcast  Ephemeris),  RSJ 
Magnitude  of  the  Difference  Vector 


Day  of 
Year  1991 

Segment 

AE(m) 

AN(m) 

AU(m) 

AH(m) 

RSS(n 

197 

A 

-23.4+  1.8 

-3.6+  1.1 

20.5  ±  1.0 

23.7+  1.9 

31.3 

B 

-1.8  ±.7 

2.9 +  .8 

12.0  ±.8 

3.4+1.1 

12.5 

198 

C 

8.8  ±  2.0 

-2.7  ±1.0 

0.7  ±.8 

9.2  ±2.2 

9.3 

D 

-5.4  ±.8 

2.3  +  .9 

17.8+  1.3 

5.9+1.1 

18.7 

E 

-7.9+1.0 

-0.3  ±  .8 

1.6  ±.8 

7.9  ±1.0 

8.0 
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V.     CONCLUSIONS 

GASP  estimated  point  positions,  produced  with  data  collected  from  three  GPS  re- 
ceivers over  seven  days,  are  compared  to  independent  Transit  Doppler  estimated  posi- 
tions. The  differences  between  the  positions  provide  a  measure  of  accuracy  for  the 
computed  GPS  positions.  Both  broadcast  and  precise  ephemerides  were  used  to  estimate 
the  GASP  point  positions.  The  precision  for  a  group  of  estimated  positions  is  deter- 
mined by  calulating  the  mean  of  the  differences  and  standard  deviations  on  the  means 
for  the  group.  Major  conclusions  are: 

1.  GASP  point  position  estimates  generated  with  the  precise  ephemeris  agree  to 
Transit  Doppler  point  position  estimates  to  within  one  meter  on  each  axis. 

2.  All  three  GPS  receivers  used  in  this  study  are  capable  of  producing  point  positions 
with  observed  errors  on  the  RSS  of  less  than  1.5  meters,  when  the  precise  ephemeris 
and  clock  states  are  utilized. 

3.  Point  positions  estimated  with  the  precise  ephemeris  are  more  accurate  and  precise 
than  positions  estimated  with  the  broadcast  ephemeris. 

4.  Selective  Availability  has  had  a  dramatic  effect  on  positions. 

The  GPS  point  positions  estimated  by  GASP  using  the  precise  ephemeris  and  clock 
states  demonstrated  good  agreement  with  positions  estimated  by  the  Transit  Doppler 
positioning  system.  The  level  of  agreement  was  within  the  one  meter  noise  level  estab- 
lished for  both  systems.  In  general,  the  Transit  Doppler  positions  were  within  the  en- 
semble of  corresponding  GPS  position  estimates.  This  would  seem  to  indicate  that  many 
of  the  systematic  errors  that  effect  the  GPS  measurements  had  been  successfully  modeled 
or  removed  by  the  GASP  algorithm.  Of  course,  systematic  errors  common  to  both  po- 
sitioning systems  may  remain. 

Accuracy  and  precision  indicate  that  all  three  GPS  receivers  examined  in  this  study 
are  capable  of  producing  geodetic-quality  point  positions.  Between-receiver  compar- 
isons of  the  estimated  positions  using  both  the  precise  and  broadcast  ephemeris  reveal 
that  the  Trimble  and  Ashtech  receivers  provide  results  that  compare  favorably  to  those 
produced  by  the  TI  4100  receiver.  This  is  not  surprising  considering  that  all  three  re- 
ceivers record  the  carrier  phase  measurements  essential  to  GASP  processing  and  high- 
precision  point  positioning.  Although  the  carrier  phase  observable  recorded  by  the  three 
receivers  is  the  same,  the  pseudorange  observable  exhibits  some  important  differences. 
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Pseudoranges  are  observed  from  the  LI  signal  using  the  C/A-code  for  the  Ashtech 
and  Trimble  receivers.  The  TI  4100  receiver  observes  pseudoranges  on  both  the  LI  and 
L2  signals  using  the  P-code.  Consequently,  a  dual  frequency  ionospheric  correction  can 
not  be  computed  for  pseudoranges  observed  with  the  Ashtech  and  Trimble  receivers. 
This,  along  with  higher  noise  levels  associated  with  C/A  code  measurements,  caused  a 
high  incidence  of  data  rejection  in  preliminary  processing  runs.  Relaxing  the 
pseudorange  edit  tolerance  to  ten  meters  allowed  the  acceptance  of  a  much  higher  per- 
centage of  observations.  Even  so,  rejection  figures  for  these  receivers  were  still  higher 
than  the  TI  4100  receiver. 

Though  a  higher  percentage  of  observations  were  rejected  with  the  Ashtech  and 
Trimble  receivers,  a  high  degree  of  accuracy  and  precision  was  maintained.  The  high 
accuracy  and  precision  levels  observed  with  these  receivers  may  be  attributed  to  their 
increased  channel  capacity.  The  ability  to  collect  and  subsequently  to  process  more  si- 
multaneous satellite  carrier  phase  measurements  provides  greater  geometric  diversity. 
More  measurements  per  epoch  and  the  incorporation  of  more  satellites  into  the  proc- 
essing run  provides  increased  system  redundancy.  These  two  factors  act  to  strengthen 
the  estimated  position  solution.  Position  solutions  produced  by  the  TI  4100  receiver 
should  also  be  improved  now  that  more  satellites  are  allowed  per  processing  run.  The 
benefits  gained  from  these  additional  measurements  will  be  nullified  if  too  many  obser- 
vations are  rejected  in  the  data  editing  process.  Thus,  it  may  be  desirable  in  future  tests 
to  determine  an  optimal  pseudorange  edit  tolerance. 

Point  positions  estimated  using  precise  ephemerides  and  clock  states  are  superior  to 
the  solutions  produced  with  the  broadcast  ephemeris  and  clock  states.  The  results  ob- 
tained with  the  three  receivers  at  the  Beach  Lab  demonstrate  that,  overall,  better  posi- 
tions estimates  were  provided  with  the  precise  ephemeris  and  clock  states.  This  is  clearly 
the  case  at  the  LOBOS3  site  where  significantly  better  results  were  achieved  with  the 
precise  ephemeris.  These  results,  together  with  results  produced  in  past  studies,  demon- 
strate that  positions  produced  with  the  precise  ephemeris  are  generally  more  reliable, 
accurate,  and  consistent  than  those  produced  with  the  broadcast  ephemeris  [Ref.  2:  p. 
497]. 

The  recent  reactivation  of  Selective  Availabilty  has  accomplished  its  intended  pur- 
pose. For  the  results  reported  in  this  study,  SA  has  severely  degraded  the  accuracy  and 
repeatability  of  point  position  estimates  that  were  produced  with  the  broadcast 
ephemeris.  The  effects  of  SA  on  position  estimates  using  the  precise  ephemeris  were  not 
assessed  in  this  investigation.  This  may  be  a  topic  for  future  consideration. 
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The  RINEX  modifications  to  the  GASP  program  now  allow  it  to  utilize  data  col- 
lected with  a  greater  variety  of  GPS  receivers.  This  will  permit  point  position  determi- 
nations from  a  multitude  of  additional  sources  and  will  enable  DMA  to  augment  its 
point  positioning  data  base.  Only  data  that  was  converted  to  the  original  RINEX  format 
was  used  in  these  tests.  Additional  tests  should  be  conducted  with  data  presented  in  the 
RINEX  2  format  to  ensure  that  the  program  code  designed  to  process  this  data  is  free 
of  error. 

Additional  tests  should  be  conducted  to  verify  the  validity  of  the  conclusions  pre- 
sented here.  Tests,  similiar  in  design  but  extending  over  longer  time  periods  (weeks  as 
opposed  to  days),  would  provide  a  more  satisfactory  indication  of  the  long-term  repeat- 
ability for  the  different  receivers. 
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APPENDIX  A.     KEPLERIAN  ELEMENTS 

Six  Keplerian  elements  and  a  reference  time  are  necessary  to  completely  describe  the 
satellite  orbit  and  the  position  of  the  satellite  in  the  orbit.  The  six  Keplerian  elements 
given  in  the  broadcast  ephemeris  are: 

•  Right  ascension  of  the  ascending  node  (Q.) 

•  Inclination  (i) 

•  Argument  of  perigee  (a>) 

•  Semi-major  axis  of  the  elliptical  orbit  (a) 

•  Eccentricity  of  the  orbit  (e) 

•  An  element  describing  the  position  of  the  satellite  on  the  orbital  ellipse  (e.g.,  Mean 
anomaly  (M))  which  is  a  function  of  time 

The  ephemeris  reference  time  is  designated  (/„).  The  five  Keplerian  elements  that  de- 
scribe the  orbit  are  Q,  i,  co,  a,  and  e.  M  gives  the  position  of  the  satellite  in  its  orbit  at 
a  time  t.  See  Figure  33  for  a  graphic  representation. 
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Other  Broadcast  Ephemeris  parameters  that  are  transmitted  in  the  Navigation  Message 
describe  the  deviations  of  the  satellite  motion  from  the  smooth  ellipse  defined  by  the  six 
Keplerian  elements  refered  to  above.  These  perturbing  terms  include: 

•  Mean  motion  difference  (An) 

•  Rate  of  right  ascension  (Q) 

•  Rate  of  inclination  (/) 
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•  Corrections  to  argument  of  latitude  (CueCut) 

•  Corrections  to  orbital  radius  (CrcC„) 

•  Corrections  to  inclination  {ClcCls) 

In  order  to  transform  the  Keplerian  orbital  parameters  to  coordinates  in  the  con- 
ventional terrestrial  system,  the  following  computations  must  be  performed: 

•  find  time  since  reference  epoch 

tk  =  t  —  toe 

•  solve  for  mean  anomaly  at  tk 


a 


where  the  gravitational  constant 
H  =  3.986005  x  1014 


sec 


• 


solve  Kepler's  equation  iteratively  for  eccentric  anomaly,  Ek 

Mk  =  Ek  —  e  sin  Ek 

•  compute  true  anomaly, 

,    ,  -if  r.     r    ( sin  Ek)    1 

•  compute  the  argument  of  latitude, 

uk  =  w  +  fk  +  Cuc  cos  2(<u  +fk)  +  Cus  sin  2(<"  +  fk) 

•  compute  the  orbital  radius, 

rk  =  a  (1  —  e  cos  E^  +  Cre  cos  2(cw  +.4)  -1-  C„  sin  2(<w  -l-,^) 

compute  the  orbital  inclination, 

k  —  'o  +  "a  +  Qe cos  2(<w  +./*)  +  Qi  s"1 2(w  +A) 
compute  the  longitude  of  the  ascending  node, 

^  =  Q0  +  (Q  -  coe)  tk  -  a)e  toe 

where  the  mean  earth  rotation  rate  is 
w,  =  7.2921 15147  xl0-s -If-- 
then compute  the  CT  (earth-fixed)  coordinates, 

Xk  =  rk  cos  uk  cos  Xk  —  rk  sin  uk  cos  ik  sin  A* 

yA  =  rk  cos  wA  sin  kk  +  rk  sin  wA  cos  ik  cos  A^ 
ZA  =  rk  sin  wA  sin  ik 
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APPENDIX  B.     LEAST  SQUARES  IN  GPS 

The  contents  of  this  appendix  follow  the  development  in  Uotila  [Ref.  16].  This  de- 
velopment focuses  on  the  general,  nonlinear,  observation  equation  model.  This  is  the 
model  utilized  by  GASP  in  performing  its  batch  least  squares  adjustment. 

A.  THE  MATHEMATICAL  MODEL 

Using  matrix  notation  the  mathematical  model  is  given  by, 

La  =  F{X°) 

La  =  FiX0) 

Lb  -  e  =  F{X°) 

where,  L"  =  Theoretical  values  of  the  observed  quantities 
X=  Theoretical  values  of  the  paramters 

A 

La  =  Estimates  or  adjusted  values  of  the  observed  quantities 

A 

X—  Estimates  or  adjusted  values  of  the  parameters 
Lh  =  Observed  values  of  the  observables 
e  =  "True"  errors 

B.  OBSERVATION  EQUATIONS 

Prior  to  employing  the  least  squares  method,  a  nonlinear  function  is  linearized  by 
expanding  in  a  Taylor  series  and  eliminating  terms  second  order  and  higher.  For  a  set 
of  nonlinear  functions  in  matrix  form  we  get, 


dX   ■*•-*"  V*        A  'T       -'VUT  d][a 


Lb-s  =  F(Xa)  =  F(X°)  +  -^\xO=:x0(r-X0)  +  -  =  F(X0)  +  ^\x.  =  x0X  + 

dF 


a=!xt1x°=x° 

where,  X0  are  approximate  values  for  the  parameters  and  is  the  point  about  which  ex- 
pansion is  done,  and 

A   is  called  the  design  matrix  and  it  has  the  form 
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A  = 


dx] 

dx2 

3/i 

dxu 

Vi 

to 

¥2 

8Xy 

dx2 

dxu 

dfn 

k 

k 

dxx       dx2 


dx„ 


where,  n  is  the  number  of  observations,  and 

u  is  the  number  of  equation  parameters.    We  will  assume  that  redundancy  exists, 
i.e.,  the  number  of  observations  is  greater  than  the  number  of  unknown  parameters. 
Using  the  relations: 

e  =  AX+F(X°)-Lb 

V=AX+F(X°)-Lb 

where,  V  is  a  vector  of  differences  between  the  observed  values  and  the  estimated  values 
(V  is  refered  to  as  the  residual  vector)  and 

X=X°-X° 


and  from, 


X-fr-X0 


F{X°)  =  L° 


L°-Lb  =  L 


where,  L°  are  the  computed  observations  evaluated  at  X°  we  get, 


V=AX+L 


These  are  the  observation  equations. 


C.    MINIMUM  VARIANCE  SOLUTION 

To  find  a  unique  set  of  parameter  estimates  some  condition  must  be  imposed  on  the 
residuals.  By  minimizing  the  sum  of  the  squares  of  the  weighted  residuals,  we  obtain  a 
best  estimate  for  the  parameters  in  a  least  squares  sense. 
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The  weight  matrix  denoted  P  is  given  as, 


P  =  °l*lX 


where,  L^1  is  the  inverse  of  the  observation  variance-covariance  matrix  and 
cr0  is  the  a  priori  variance  of  unit  weight. 
The  weight  matrix  provides  a  way  to  assign  relative  importance  to  various  observa- 
tions. The  observation  variance-covariance  matrix  lets  us  express  random  errors  or 
noise  associated  with  the  observations  as  variances  and  relationships  between  observa- 
tions as  co variances.  It  is  through  the  variance-covariance  matrix  that  knowledge  about 
the  observations  is  propagated.  This  allows  the  variances  and  covariances  for  the  esti- 
mated parameters  to  be  determined  (i.e.,  it  gives  an  estimate  of  the  noise  in  the  solution). 
The  variances  make  up  the  diagonal  elements  of  the  matrix  and  the  off-diagonal  ele- 
ments are  the  covariances.  The  variance-covariance  matrix  may  be  written, 


aL\     aL\  LI 
2 
aL2  L\     °L2 


°Ln  L\  aLn  11 


°L\lM 
aL2Ln 


2 


where,  au  =  aJI  i.e.,  a  symmetric  matrix.    The  variance-covariance  matrix  is  diagonal  if 
the  observations  are  uncorrelated. 

For  the  minimum  variance  solution,  we  want  to  minimize  V7  P  V  which  is  the  sum 
of  the  squares  of  the  weighted  residuals  in  matrix  notation.  By  some  algebraic  manipu- 
lation we  can  show, 


V1  PV^XA1  PAX+2L1  PAX+LTPL 

A 

Minimizing  by  partially  differentiating  with  respect  to  X  gives, 


1      d(VTPV) 


A 

dX 


=  AJ  PAX+AJ  PL  =  0 


This  is  known  as  the  normal  equation.  Now  solving  for  X  we  get, 
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X=-{AT PA)~^  AT PL 

A 

where,  X  is  the  adjustment  vector  that  is  applied  to  X*  giving  an  improved  estimate  for 
the  parameters.  A  number  of  iterations  may  be  necessary  before  X  reaches  some  ac- 
ceptable level,  that  is  it  sufficiently  converges. 

D.     A  POSTERIORI  VARIANCE  OF  UNIT  WEIGHT 

The  sum  of  the  squares  of  the  residuals  V  P  V  are  found  by, 

VTPV=LTPL  +  XrATPL 
The  a  posteriori  variance  of  unit  weight  may  be  calculated, 

/V2        VT  P  V 


On  = 


where,  n  —  u  is  the  degree  of  freedom. 

E.  VARIANCE-COVARIANCE  MATRIX  OF  ADJUSTED  PARAMETERS 

The  variance-co variance  matrix  for  the  adjusted  parameters  is  given  by, 

I,fr  =  (AT  P  A)'1  AT  PI,LP  A(AT  P  A)~l 

but  with  P  =  a\  1^  this  reduces  to, 

I.^  =  ol(ATPA)-1 

F.  VARIANCE-COVARIANCE  MATRIX  OF  THE  ADJUSTED,  OBSERVED 
QUANTITIES 

This  matrix  is  given  by, 

Ill.  =  a20A{ATPA)~lAT 

G.  SUMMARY 

To  summarize  the  theoretical  aspects  of  the  least  squares  method  using  the  nonlin- 
ear, observation  equation  model, 

•  The  nonlinear  system  of  equations  is  linearized  by  expanding  in  a  Taylor  series  and 
nonlinear  terms  are  truncated  (Note:  the  number  of  equations  must  be  greater  than 
or  equal  to  the  number  of  parameters), 
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•  A  suitable  noise  model  is  chosen  as  represented  by  the  observation  variance- 
covariance  matrix  and  a  priori  variance  of  unit  weight,  and 

•  the  sum  of  the  squares  of  the  residuals  are  minimized  to  yield  a  ininimum  variance 
solution 

The  operational  procedures  to  estimate  values  for  the  unknown  parameters  are, 

•  Develop  a  mathematical  model  that  expresses  the  observations  as  a  function  of 
some  unknown  parameters, 

•  select  appropriate  approximate  values  for  the  parameters  X°  (e.g.,  the  receiver  co- 
ordinates), 

•  Form  the  vector  of  computed  observations  L°  by  evaluating  the  mathematical 
model  at  X°, 

•  Determine  L,  which  is  the  difference  vector  of  the  computed  observations  minus 
the  actual  observations, 

•  Form  the  weight  matrix  P  from  the  observation  variance-covariance  matrix  and  the 
a  priori  variance  of  unit  weight 

•  take  the  partial  derivatives  of  the  mathematical  model  with  respect  to  the  un- 
knowns and  evalute  at  X°  to  form  the  elements  of  the  design  matrix  A, 

•  solve  for  the  change  in  parameters,  X, 

•  Compute  the  parameter  variance-covariance  matrix,  a  posteriori  variance,  the  resi- 
dual vector,  etc., 

A  A 

•  apply  X  to  the  initial  approximations  and  stop  if  X  has  suiFiciently  converged, 
otherwise  repeat  the  process  using  X"  as  the  new  X°  (i.e.,  iterate) 
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APPENDIX  C.     CONSTANTS  USED  IN  GASP 

Table  21  presents  the  constants  used  in  the  GASP  program.    The  constants  are 
identified;  and  the  common  symbols,  the  units  (if  applicable),  and  the  values  are  given. 

Table  15.     CONSTANTS  USED  IN  GASP 


Constant 

Symbol 

Units 

Value 

WGS-84  Earth's  Semi-major  Axis 

a 

km 

6378.137 

WGS-84  Inverse  Earth  Flattening 

1// 

298.257223563 

Product  of  Gravitational  Constant  and 
Earth's  Mass 

GM 

kmzl  sec2 

3.986005  x  10s 

Earth's  Rotation  Rate 

(O 

rad\  sec 

7.2921151467  x  10- 

Pi 

it 

3.1415926535897932 

Speed  of  Light 

c 

kmj  sec 

299,792.458 

Satellite  Base  Frequency 

/. 

Hz 

10.23  x  10s 

LI  Signal  Frequency 

Ju 

Hz 

154/0 

L2  Signal  Frequency 

Ju 

Hz 

120/0 

LI  Signal  Wavelength 

hi 

cm 

19.029 

L2  Signal  Wavelength 

h2 

cm 

24.421 

TI  4100  Receiver  LI  Frequency  Offset 

Hz 

-6000 

TI  4100  Receiver  L2  Frequency  Offset 

Hz 

7600 
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APPENDIX  D.     TABLES  OF  RESULTS  IN  CT  COORDINATES 

A.     BEACH  LAB  RESULTS 

The  results  previously  expressed  in  terms  of  local  east,  north,  and  up  coordinates 
are  presented  here  as  component  differences  in  conventional  terrestrial  coordinates  (X, 
Y,  and  Z)  for  the  reader  more  comfortable  with  this  reference  frame.  Tables  15  through 
17  display  these  results.  No  corresponding  target  plots  are  shown. 


Table  16.  BEACH  LAB  COMPONENT  DIFFERENCES  (AX,  AY,  AZ)  USING  BROADCAST 
EPHEMERIS:  Solution  Differences  Between  GPS  and  Transit  Doppler  Methods, 
Second  Term  is  GASP  Estimated  Formal  Error 


Day  of 
Year  1990 

Receiver 

AX(m) 

AY(m) 

AZ(m) 

338 

TI  4100 

-1.1  +  .6 

-0.6  ±  .6 

-1.0  ±.4 

Ashtech 

1.1  ±.7 

1.1 +  .4 

-1.0  ±.4 

Trimble 

1.0 +  .6 

2.6 +  .4 

-0.6  ±  .4 

339 

TI  4100 

-1.5  ±.6 

-1.7  ±.7 

-0.2  ±  .4 

Ashtech 

4.6  ±.7 

1.1 +  .5 

-1.3  ±.4 

Trimble 

4.5  ±  .7 

1.8  ±.4 

-1.4  ±.4 

340 

TI  4100 

-0.9  ±  .6 

-0.4  ±  .6 

-0.2  ±  .4 

Ashtech 

1.4  ±.6 

0.8  ±  .4 

-1.2  ±.4 

Trimble 

1.2  ±.6 

1.2  ±.4 

-0.3  ±  .4 

341 

TI  4100 

1.2  ±.5 

-0.9  ±  .5 

-0.9  ±  .4 

Ashtech 

-1.2  ±.6 

1.5  ±.4 

-0.0  +  .4 

Trimble 

0.2  ±.6 

-0.4  ±  .4 

0.1  ±.4 
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Table  17. 


BEACH  LAB  COMPONENT  DIFFERENCES  (AX,  AY,  AZ)  USING  PREC1 
EPHEMERIS:  Solution  Differences  Between  GPS  and  Transit  Doppler  Methc 
Second  Term  is  GASP  Estimated  Formal  Error 


Day  of 
Year  1990 

Receiver 

AX(m) 

AY(m) 

AZ(m) 

338 

TI  4100 

2.2 +  .6 

1.1  ±.6 

0.1 +  .4 

Ashtech 

-0.9  ±  .6 

0.9  ±  .4 

-0.9  ±  .4 

Trimble 

-1.2 +  .6 

1.7  ±.4 

-0.2  ±  .4 

339 

TI  4100 

-0.2  ±  .6 

0.6  ±.6 

-0.4  ±  .4 

Ashtech 

-0.5  ±  .7 

-0.1 +  .5 

-0.5  ±  .4 

Trimble 

0.6  ±.6 

0.9  ±  .4 

-0.9  ±  .4 

340 

TI  4100 

-0.7  ±  .6 

1.9  ±.6 

-1.0 +.4 

Ashtech 

0.8 +  .6 

1.4  ±.4 

-1.7  ±.4 

Trimble 

0.3  ±  .6 

1.2 +  .4 

-0.6  ±  .4 

341 

TI  4100 

1.3  ±.6 

-1.8  ±.6 

1.6  ±.4 

Ashtech 

0.1 +  .6 

0.5  ±  .5 

0.5  ±  .4 

Trimble 

1.0  ±.6 

-1.0 +  .4 

0.5  ±  .4 
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Table  18.     BEACH  LAB  REPEATABILITY  (AX,  AY,  AZ):  AVERAGE  OVER  DAYS  BY  RE- 
CEIVER 


Ephemeris 

Receiver 

Mean  AX  (m) 

Mean  AY  (m) 

Mean  AZ  (m) 

Broadcast 

TI  4100 

-0.6  ±  1.2 

-0.9  +  .6 

-0.6  ±  .4 

Ashtech 

1.5  +  2.4 

1.1  +  .3 

-0.9  ±  .6 

Trimble 

1.7  ±  1.9 

1.3  ±  1.3 

-0.5  ±.6 

Precise 

TI  4100 

0.7+  1.3 

0.4  ±  1.6 

0.1  ±  1.1 

Ashtech 

-0.1 +  .7 

0.7 +  .6 

-0.7 +  .9 

Trimble 

0.2  ±.9 

0.7+1.2 

-0.3  ±  .6 
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B.     LOBOS3  RESULTS 

Tables  18  and  19  display  the  LOBOS3  results  in  conventional  terrestrial  X,  Y,  and 
Z  coordinates.  Associated  target  plots  are  not  presented. 


Table  19.     LOBOS3  COMPONENT  DIFFERENCES  (AX,  AY,  AZ):     Solution  Differences 
tween  GPS  and  Transit  Doppler  Methods  for  Broadcast  and  Precise  Ephemerides,  S 
ond  Term  is  GASP  Estimated  Formal  Error 


Ephemeris 

Day  of  Year  1991 

AX(m) 

AY(m) 

AZ(m) 

Broadcast 

36 

-0.1  ±.7 

-2.7  ±  .5 

-0.2  ±  .4 

37 

-0.6  ±  .9 

-7.5  ±  .6 

-0.8  ±  .5 

38 

-2.4  ±.8 

-6.4  ±  .5 

0.8  ±  .4 

Precise 

36 

-0.7  ±  .7 

-0.5  ±  .5 

0.0  ±  .4 

37 

0.5  ±  .7 

0.1 +  .4 

-0.5  ±  .3 

38 

1.0  ±.7 

-0.8  ±  .5 

0.2 +  .4 
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Table  20.     LOBOS3  REPEATABILITY  (AX,  AY,  AZ):  AVERAGE 
OVER  DAYS 


Ephemeris 

Mean  AX  (m) 

Mean  AY  (m) 

Mean  AZ  (m) 

Broadcast 

-1.0  ±  1.2 

-5.5  ±2.5 

-0.1  +  .8 

Precise 

0.3 +  .9 

-0.4  ±  .5 

-0.1  ±.4 
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C.     BLDG  224.3  RESULTS 

Tat      20  presents  the  solution  differences  for  BLDG  224.3  in  the  CT  coordinate 
system. 


Table  21.     BLDG  224.3  COMPONENT  DIFFERENCES  (AX,  AY,  AZ):     Selective  Availab 
Test  -  GPS  Solution  Differences  Between  Previously  Determined  Coordinates  and 
ordinates  Determined  After  Activation  of  S/A  (Using  Broadcast  Ephemeris) 


Day  of 
Year  1991 

Segment 

AX(m) 

AY(m) 

AZ(m) 

197 

A 

-29.7  ±1.9 

-4.5  ±1.1 

5.8  ±.6 

B 

-5.7 +  .8 

-5.8  ±.9 

9.5  ±.5 

198 

C 

6.4  ±2.1 

-6.5  ±  .9 

-1.7  ±.7 

D 

-11.4+1.0 

-8.1  ±1.3 

12.4  ±.7 

E 

-7.5  ±1.0 

2.9  ±1.1 

0.7  ±  .5 
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APPENDIX  E.     SUMMARIES  OF  POSITIONING  RESULTS 
A.     SUMMARY  OF  TRANSIT  POSITION  RESULTS 


Station  name: 

DOP5 

Method: 

Transit  Doppler  Observations  (@DOP3) 

and  Conventional  Terrestrial  Surveying 

Period  of  Occupation: 

January  26-30,  1991 

Transit  Receiver: 

Magnavox  MX- 1502 

Satellites  Observed: 

77,105,115,124,126,128 

Passes  Accepted: 

57 

RMS  of  Residuals: 

8.0  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707273.74  -4353292.80  3781989.36 

Geodetic  Coordinates  (monument): 

(D  =  36°36'06.2347    A  = 

238°07'22.404    h  =  -23.55  meters 

Standard  Deviations  (ct>,A,h)  (meters): 

1.0   1.0   1.0  --   1  Sigma 
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Station  name: 

DOP1 

Method: 

Transit  Doppler  Observations  (@DOP3) 
and  Conventional  Terrestrial  Surveying 

Period  of  Occupation: 

January  26-30,  1991 

Transit  Receiver: 

Magna  vox  MX- 1502 

Satellites  Observed: 

77,105,115,124,126,128 

Passes  Accepted: 

57 

RMS  of  Residuals: 

8.0  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707255.23  -4353301.05  3781993.34 

Geodetic  Coordinates  (monument): 

<D  =  36°36'06.3917 

A  = 

238°07'23.211     h  =  -23.40  meters 

Standard  Deviations  (<D,A,h)  (meters): 

1.0   1.0   1.0  -   1  Sigma 

Station  name: 

DOP2 

Method: 

Transit  Doppler  Observations  (@DOP3) 
and  Conventional  Terrestrial  Surveying 

Period  of  Occupation: 

January  26-30,  1991 

Transit  Receiver: 

Magna  vox  MX- 1502 

Satellites  Observed: 

77,105,115,124,126,128 

Passes  Accepted: 

57 

RMS  of  Residuals: 

8.0  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707259.43  -4353297.54  3781994.29 

Geodetic  Coordinates  (monument): 

<D  =  36°36'06.4311 

A  = 

238°07'22.9936    h  =  -23.45  meters 

Standard  Deviations  (0,A,h)  (meters) 

1.0   1.0   1.0  --   1  Sigma 
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Station  name: 

LOBOS3 

Method: 

Transit  Doppler  Observations 

Period  of  Occupation: 

January  31- 

February  4,  1991 

Transit  Receiver: 

Magnavox  MX- 1502 

Satellites  Observed: 

77,105,115,124,126,128 

Passes  Accepted: 

62 

RMS  of  Residuals: 

9.0  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2693399.38 

-4346211.98  3799864.81 

Geodetic  Coordinates  (monument): 

<D  =  36°48'09.683 

A  = 

238°12'46.841 

h  =  —30.45  meters 

Standard  Deviations  (<D,A,h)  (meters) 

1.0   1.0   1.0 

-   1  Sigma 
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B.    SUMMARY  OF  GPS  POSITION  RESULTS 
1.     Beach  Lab  Sites 


Station  name: 

DOP5 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0J 

Ephemerides: 

Broadcast 

Satellite  Clock  States: 

Broadcast 

Date  of  Occupation: 

Day  338,  1990 

GPS  Receiver: 

TI  4100  BEPP/CORE  (versions  3.7/5.1) 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6.2  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,13,18,19 

Final  Number  of  GASP  Observables: 

1030 

Percentage  of  Data  Rejected: 

2.2  % 

RMS  of  Residuals: 

4.709  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707274.8418  -4353293.3659  3781988.4006 

Standard  Deviations  (X,Y,Z)  (meters) 

: 

0.6310  0.5967  0.4054 

Geodetic  Coordinates  (monument): 

<D  =  36°36'06.1889    A  =  : 

238°07'22.3782    h  =  -23.268  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.5403  meters 
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Station  name: 

DOP1 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Broadcast 

Satellite  Clock  States: 

Broadcast 

Date  of  Occupation: 

Day  339,  1990 

GPS  Receiver: 

TI  4100  BEPP/CORE  (versions  3.7/5.1) 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6.2  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,13,18,19 

Final  Number  of  GASP  Observables: 

1017 

Percentage  of  Data  Rejected: 

1.2  % 

RMS  of  Residuals: 

4.794  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707256.7432  -4353302.7647  3781993.1333 

Standard  Deviations  (X,Y,Z)  (meters) 

0.6480  0.6657  0.4016 

Geodetic  Coordinates  (monument): 

<D  =  36°36'06.3426    A  =  : 

tt8°07'23. 1963    h  =  -21 .712  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.3673  meters 

HI 


Station  name: 

DOP5 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0; 

Ephemerides: 

Broadcast 

Satellite  Clock  States: 

Broadcast 

Date  of  Occupation: 

Day  340,  1990 

GPS  Receiver: 

TI  4100  BEPP/CORE  (versions  3.7/5.1) 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

7.1  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,13,18,19 

Final  Number  of  GASP  Observables: 

1185 

Percentage  of  Data  Rejected: 

2.3  % 

RMS  of  Residuals: 

4.584  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707276.6237  -4353293.1626  3781989.1370 

Standard  Deviations  (X,Y,Z)  (meters) 

0.5496  0.5643  0.3848 

Geodetic  Coordinates  (monument): 

<D  =  36°36'06.1932    A  =  : 

Z38°07'22.3130    h  =  -22.212  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.5873  meters 
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Station  name: 

DOP5 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Broadcast 

Satellite  Clock  States: 

Broadcast 

Date  of  Occupation: 

Day  341,  1990 

GPS  Receiver: 

TI  4100  BEPP/CORE  (versions  3.7/5.1) 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6.5  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,13,18,19 

Final  Number  of  GASP  Observables: 

690 

Percentage  of  Data  Rejected: 

38.2  % 

RMS  of  Residuals: 

3.306  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707272.5806  -4353293.6093  3781988.4578 

Standard  Deviations  (X,Y,Z)  (meters) 

0.5463  0.5307  0.3803 

Geodetic  Coordinates  (monument): 

<D  =  36°36'06.2095    A  =  238°07'22.4606    h  -  -24.027  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.5103  meters 
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Station  name: 

DOP1 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Broadcast 

Satellite  Clock  States: 

Broadcast 

Date  of  Occupation: 

Day  338,  1990 

GPS  Receiver: 

Ashtech  LD  XII 

Receiver  software  version 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

5.6  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,15,16,18,19 

Final  Number  of  GASP  Observables: 

1733 

Percentage  of  Data  Rejected: 

9.3  % 

RMS  of  Residuals: 

3.336  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707254.1323  -4353299.9074  3781992.3652 

Standard  Deviations  (X,Y,Z)  (meters) 

0.6568  0.4494  0.4107 

Geodetic  Coordinates  (monument): 

<D  =  36°36'06.3962    A  =  238°07'23.2248    h  =  -25.225  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.3313  meters 
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Station  name: 

DOP5 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0] 

Ephemerides: 

Broadcast 

Satellite  Clock  States: 

Broadcast 

Date  of  Occupation: 

Day  339,  1990 

GPS  Receiver: 

Ashtech  LD  XII 

Receiver  software  version 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,15,16,18,19 

Final  Number  of  GASP  Observables: 

1079 

Percentage  of  Data  Rejected: 

19.8% 

RMS  of  Residuals: 

4.590  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707269.1762  -4353291.6786  3781988.0816 

Standard  Deviations  (X,Y,Z)  (meters) 

0.7447  0.4817  0.4241 

Geodetic  Coordinates  (monument): 

O  =  36°36'06.2662    A  =  238°07'22.5359    h  =  -27.01 1  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.3373  meters 
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Station  name: 

DOP1 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0] 

Ephemerides: 

Broadcast 

Satellite  Clock  States: 

Broadcast 

Date  of  Occupation: 

Day  340,  1990 

GPS  Receiver: 

Ashtech  LD  XII 

Receiver  software  version 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6.8  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,15,16,18,19 

Final  Number  of  GASP  Observables:     1297 

Percentage  of  Data  Rejected: 

9.7  % 

RMS  of  Residuals: 

4.481  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707253.7859  -4353300.2267  3781992.1104 

Standard  Deviations  (X,Y,Z)  (meters): 

0.6444  0.4487  0.4048 

Geodetic  Coordinates  (monument): 

<D  =  36°36'06.3879    A  = 

=  238°07'23.2434    h  =  -25.306  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.2483  meter s 
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Station  name: 

DOP1 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Broadcast 

Satellite  Clock  States: 

Broadcast 

Date  of  Occupation: 

Day  341,  1990 

GPS  Receiver: 

Ashtech  LD  XII 

Receiver  software  version 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6.9  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,15,16,18,19 

Final  Number  of  GASP  Observables: 

1255 

Percentage  of  Data  Rejected: 

13.6% 

RMS  of  Residuals: 

4.259  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707256.4461  -4353299.5540  3781993.3044 

Standard  Deviations  (X,Y,Z)  (meters) 

0.6027  0.4345  0.3941 

Geodetic  Coordinates  (monument): 

O=36°36'06.4029    A  =  238°07'23.1382    h  -  -23.925  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.4823  meters 
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Station  name: 

DOP2 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Broadcast 

Satellite  Clock  States: 

Broadcast 

Date  of  Occupation: 

Day  338,  1990 

GPS  Receiver: 

Trimble  4000ST 

Receiver  software  version  4.3X 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,15,16,18,19 

Final  Number  of  GASP  Observables:     1151 

Percentage  of  Data  Rejected: 

11.1% 

RMS  of  Residuals: 

4.014  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707258.5156  -4353294.9312  3781993.6745 

Standard  Deviations  (X,Y,Z)  (meters): 

0.6452  0.4227  0.3927 

Geodetic  Coordinates  (monument): 

■•' 

<D  =  36°36'06.4673    A  = 

=  238°07'22.9693     h  =  -25.978  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.3313  meters 
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Station  name: 

DOP2 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Broadcast 

Satellite  Clock  States: 

Broadcast 

Date  of  Occupation: 

Day  339,  1990 

GPS  Receiver: 

Trimble  4000ST 

Receiver  software  version  4.3X 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,15,16,18,19 

Final  Number  of  GASP  Observables: 

1073 

Percentage  of  Data  Rejected: 

18.0% 

RMS  of  Residuals: 

3.929  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707254.9181  -4353295.7659  3781992.9337 

Standard  Deviations  (X,Y,Z)  (meters) 

0.6625  0.4241  0.3944 

Geodetic  Coordinates  (monument): 

<D  =  36°36'06.4710    A  =  238°07'23.1099    h  -  -27.376  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.3023  meters 
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Station  name: 

DOP2 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Broadcast 

Satellite  Clock  States: 

Broadcast 

Date  of  Occupation: 

Day  340,  1990 

GPS  Receiver: 

Trimble  4000ST 

Receiver  software  version  4.3X 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

7.5  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,15,16,18,19 

Final  Number  of  GASP  Observables:     1 282 

Percentage  of  Data  Rejected: 

18.8% 

RMS  of  Residuals: 

3.896  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707258.2384  -4353296.3532  3781994.0224 

Standard  Deviations  (X,Y,Z)  (meters): 

0.5624  0.4229  0.3724 

Geodetic  Coordinates  (monument): 

<D  =  36°36'06.4558    A  = 

=  238°07'23.0090    h  =  -24.919  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.2743  meters 
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Station  name: 

DOP2 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Broadcast 

Satellite  Clock  States: 

Broadcast 

Date  of  Occupation: 

Day  341,  1990 

GPS  Receiver: 

Trimble  4000ST 

Receiver  software  version  4.3X 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6.7  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,15,16,18,19 

Final  Number  of  GASP  Observables: 

1280 

Percentage  of  Data  Rejected: 

13.5% 

RMS  of  Residuals: 

4.096  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707259.2464  -4353297.8957  3781994.4146 

Standard  Deviations  (X,Y,Z)  (meters) 

0.6003  0.4309  0.3825 

Geodetic  Coordinates  (monument): 

O  =  36°36'06.4304    A  =  238°07'23.0073     h  =  -23.206  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.3933  meters 
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Station  name: 

DOP5 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Precise  WGS  84  EF##90336 

Satellite  Clock  States: 

Precise  PC90336 

Date  of  Occupation: 

Day  338,  1990 

GPS  Receiver: 

TI  4100  BEPP/CORE  (versions  3.7/5.1) 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6.2  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,13,18,19 

Final  Number  of  GASP  Observables: 

1033 

Percentage  of  Data  Rejected: 

1.9  % 

RMS  of  Residuals: 

4.717  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707271.5459  -4353291.7286  3781989.5013 

Standard  Deviations  (X,Y,Z)  (meters) 

i 

0.6306  0.5974  0.4060 

Geodetic  Coordinates  (monument): 

<D  =  36°36'06.2781     A  =  \ 

238°07'22.4560    h  =  -25.125  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.5403  meters 
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Station  name: 

DOP1 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0] 

Ephemerides: 

Precise  WGS84  EF##90336 

Satellite  Clock  States: 

Precise  PC90336 

Date  of  Occupation: 

Day  339,  1990 

GPS  Receiver: 

TI  4100  BEPP/CORE  (versions  3.7/5.1) 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6.2  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,13,18,19 

Final  Number  of  GASP  Observables: 

1021 

Percentage  of  Data  Rejected: 

0.8  % 

RMS  of  Residuals: 

4.361  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707255.4509  -4353300.4770  3781992.9055 

Standard  Deviations  (X,Y,Z)  (meters) 

: 

0.5886  0.6084  0.3707 

Geodetic  Coordinates  (monument): 

O  -  36°36'06.3875    A  =  : 

238°07'23.1918    h  =  -23.955  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.3673  meters 
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Station  name: 

DOP5 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0] 

Ephemerides: 

Precise  WGS  84  EF##90336 

Satellite  Clock  States: 

Precise  PC90336 

Date  of  Occupation: 

Day  340,  1990 

GPS  Receiver: 

TI  4100  BEPP/CORE  (versions  3.7/5.1) 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

7.1  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,13,18,19 

Final  Number  of  GASP  Observables: 

1182 

Percentage  of  Data  Rejected: 

2.6  % 

RMS  of  Residuals: 

4.563  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707274.3957  -4353290.8606  3781988.3570 

Standard  Deviations  (X,Y,Z)  (meters) 

0.5548  0.5665  0.3870 

Geodetic  Coordinates  (monument): 

<D  =  36°36'06.2334    A  =  238°07'22.3402    h  =  -25.191  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.5873  meters 
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Station  name: 

DOP5 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Precise  WGS  84  EF##90336 

Satellite  Clock  States: 

Precise  PC90336 

Date  of  Occupation: 

Day  341,  1990 

GPS  Receiver: 

TI  4100  BEPP/CORE  (versions  3.7/5.1) 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6.5  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,13,18,19 

Final  Number  of  GASP  Observables: 

1090 

Percentage  of  Data  Rejected: 

2.4  % 

RMS  of  Residuals: 

4.506  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707272.3968  -4353294.6322  3781990.9281 

Standard  Deviations  (X,Y,Z)  (meters) 

0.5813  0.5608  0.3847 

Geodetic  Coordinates  (monument): 

<D  =  36°36'06.2589    A  =  : 

>38°07'22.4886    h  =  -21.934  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.5103  meters 
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Station  name: 

DOP1 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0] 

Ephemerides: 

Precise  WGS  84  EF##90336 

Satellite  Clock  States: 

Precise  PC90336 

Date  of  Occupation: 

Day  338,  1990 

GPS  Receiver: 

Ashtech  LD  XII 

Receiver  software  version 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

5.6  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,15,16,18,19 

Final  Number  of  GASP  Observables: 

1737 

Percentage  of  Data  Rejected: 

9.1  % 

RMS  of  Residuals: 

3.265  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707256.1438  -4353300.1852  3781992.4439 

Standard  Deviations  (X,Y,Z)  (meters) 

'- 

0.6395  0.4393  0.4035 

Geodetic  Coordinates  (monument): 

0>  =  36°36'06.3732    A  =  : 

238°07'23.1619     h  =  -24. 136  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.3313  meters 
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Station  name: 

DOP5 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0; 

Ephemerides: 

Precise  WGS  84  EF##90336 

Satellite  Clock  States: 

Precies  PC90336 

Date  of  Occupation: 

Day  339,  1990 

GPS  Receiver: 

Ashtech  LD  XII 

Receiver  software  version 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,15,16,18,19 

Final  Number  of  GASP  Observables: 

1092 

Percentage  of  Data  Rejected: 

18.8% 

RMS  of  Residuals: 

4.347  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707274.2183  -4353292.9093  3781988.8495 

Standard  Deviations  (X,Y,Z)  (meters) 

0.7005  0.4521  0.4073 

Geodetic  Coordinates  (monument): 

<D  =  36°36'06.2144    A  =  238°07'22.3898    h  =  -23.576  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.3373  meters 
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Station  name: 

DOP1 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0J 

Ephemerides: 

Precise  WGS  84  EF##90336 

Satellite  Clock  States: 

Precise  PC90336 

Date  of  Occupation: 

Day  340,  1990 

GPS  Receiver: 

Ashtech  LD  XII 

Receiver  software  version 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6.8  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,15,16,18,19 

Final  Number  of  GASP  Observables: 

1297 

Percentage  of  Data  Rejected: 

9.7  % 

RMS  of  Residuals: 

4.338  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707254.4187  -4353299.6465  3781991.6860 

Standard  Deviations  (X,Y,Z)  (meters) 

0.6276  0.4352  0.3971 

Geodetic  Coordinates  (monument): 

(D  =  36°36'06.3799    A  =  238°07'23.2094    h  =  -25.686  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.2483  meters 
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Station  name: 

DOP1 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0J 

Ephemerides: 

Precise  WGS  84  EF##90336 

Satellite  Clock  States: 

Precise  PC90336 

Date  of  Occupation: 

Day  341,  1990 

GPS  Receiver: 

Ashtech  LD  XII 

Receiver  software  version 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6.9  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,15,16,18,19 

Final  Number  of  GASP  Observable! 

>:     1263 

Percentage  of  Data  Rejected: 

13.1% 

RMS  of  Residuals: 

4.420  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707255.1764  -4353300.5816  3781993.8009 

Standard  Deviations  (X,Y,Z)  (meters): 

0.6308  0.4529  0.4070 

Geodetic  Coordinates  (monument): 

O  =  36o36'06.4119    A  = 

=  238°07'23.2034    h  =  -23.466  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.4823  meters 
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Station  name: 

DOP2 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0] 

Ephemerides: 

Precise  WGS  84  EF##90336 

Satellite  Clock  States: 

Precise  PC90336 

Date  of  Occupation: 

Day  338,  1990 

GPS  Receiver: 

Trimble  4000ST 

Receiver  software  version  4.3X 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,15,16,18,19 

Final  Number  of  GASP  Observables:     1 1 64 

Percentage  of  Data  Rejected: 

11.4% 

RMS  of  Residuals: 

3.807  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707260.5913  -4353295.8412  3781994.0927 

Standard  Deviations  (X,Y,Z)  (meters): 

0.6084  0.3998  0.3763 

Geodetic  Coordinates  (monument): 

O  =  36°36'06.4420    A  = 

=  238°07'22.9 1 77    h  =  -24.228  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.3313  meters 
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Station  name: 

DOP2 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Precise  WGS  84  EF##90336 

Satellite  Clock  States: 

Precise  PC90336 

Date  of  Occupation: 

Day  339,  1990 

GPS  Receiver: 

Trimble  4000ST 

Receiver  software  version  4.3X 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,15,16,18,19 

Final  Number  of  GASP  Observables: 

1089 

Percentage  of  Data  Rejected: 

16.8% 

RMS  of  Residuals: 

3.715  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707258.8586  -4353296.6201  3781993.4158 

Standard  Deviations  (X,Y,Z)  (meters) 

0.6186  0.3956  0.3804 

Geodetic  Coordinates  (monument): 

cp  =  36°36'06.4293    A  =  : 

238°07'22.9934    h  =  -24.836  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.3023  meters 
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Station  name: 

DOP2 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Precise  WGS  84  EF##90336 

Satellite  Clock  States: 

Precise  PC90336 

Date  of  Occupation: 

Day  340,  1990 

GPS  Receiver: 

Trimble  4000ST 

Receiver  software  version  4.3X 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

7.5  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,15,16,18,19 

Final  Number  of  GASP  Observables: 

1281 

Percentage  of  Data  Rejected: 

18.8% 

RMS  of  Residuals: 

3.793  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707259.1230  -4353296.3184  3781993.6501 

Standard  Deviations  (X,Y,Z)  (meters) 

0.5582  0.4141  0.3703 

Geodetic  Coordinates  (monument): 

■• 

<D  =  36°36'06.4377    A  =  238°07'22.9780    h  =  -24.789  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.2743  meters 
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Station  name: 

DOP2 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0J 

Ephemerides: 

Precise  WGS  84  EF##90336 

Satellite  Clock  States: 

Precise  PC90336 

Date  of  Occupation: 

Day  341,  1990 

GPS  Receiver: 

Trimble  4000ST 

Receiver  software  version  4.3X 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

6.7  Hours 

PRN  Numbers  Tracked: 

2,6,9,11,12,15,16,18,19 

Final  Number  of  GASP  Observables: 

1285 

Percentage  of  Data  Rejected: 

13.1% 

RMS  of  Residuals: 

4.214  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2707258.4303  -4353298.5061  3781994.7611 

Standard  Deviations  (X,Y,Z)  (meters) 

0.6086  0.4354  0.3854 

Geodetic  Coordinates  (monument): 

<D  =  36°36'06.4377    A  =  238°07'23.0481     h  =  -22.929  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.3933  meters 
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2.     LOBOS3  Site 


Station  name: 

LOBOS3 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0] 

Ephemerides: 

Broadcast 

Satellite  Clock  States: 

Broadcast 

Date  of  Occupation: 

Day  36,  1991 

GPS  Receiver: 

TI  4100  BEPP/CORE  (versions  3.7/5.1) 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

8.6  Hours 

PRN  Numbers  Tracked: 

3,6,11,12,13,16,17,20 

Final  Number  of  GASP  Observables: 

964 

Percentage  of  Data  Rejected: 

2.0  % 

RMS  of  Residuals: 

4.847  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2693399.4335  -4346214.6434  3799864.5746 

Standard  Deviations  (X,Y,Z)  (meters) 

0.6637  0.5214  0.3612 

Geodetic  Coordinates  (monument): 

<D  =  36°48'09.6322    A  =  : 

238°12'46.8955    h  -  -28.749  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.5673  meters 
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Station  name: 

LOBOS3 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Broadcast 

Satellite  Clock  States: 

Broadcast 

Date  of  Occupation: 

Day  37,  1991 

GPS  Receiver: 

TI  4100  BEPP/CORE  (versions  3.7/5.1) 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

10.8  Hours 

PRN  Numbers  Tracked: 

3,6,11,12,13,16,17,20 

Final  Number  of  GASP  Observable? 

:     1361 

Percentage  of  Data  Rejected: 

0.9  % 

RMS  of  Residuals: 

6.786  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2693400.0239  -4346219.4479  3799864.0505 

Standard  Deviations  (X,Y,Z)  (meters): 

0.8974  0.6011  0.4497 

Geodetic  Coordinates  (monument): 

fl>  =  36°48'09.5332    A  = 

=  238°12'46.9773    h  =  -25.544  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1 .5673  meters 
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Station  name: 

LOBOS3 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Broadcast 

Satellite  Clock  States: 

Broadcast 

Date  of  Occupation: 

Day  38,  1991 

GPS  Receiver: 

TI  4100  BEPP/CORE  (versions  3.7/5.1) 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

10.8  Hours 

PRN  Numbers  Tracked: 

3,6,11,12,13,16,17,20 

Final  Number  of  GASP  Observables: 

1293 

Percentage  of  Data  Rejected: 

0.8  % 

RMS  of  Residuals: 

5.733  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2693401.7366  -4346218.3372  3799865.5840 

Standard  Deviations  (X,Y,Z)  (meters) 

0.7575  0.5162  0.3880 

Geodetic  Coordinates  (monument): 

O  =  36°48'09.5738    A  =  238°12'46.8950    h  =  -24.659  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.5673  meters 
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Station  name: 

LOBOS3 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Precise  WGS  84  EF##91034 

Satellite  Clock  States: 

Precise  PC91034 

Date  of  Occupation: 

Day  36,  1991 

GPS  Receiver: 

TI  4100  BEPP/CORE  (versions  3.7/5.1) 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

8.6  Hours 

PRN  Numbers  Tracked: 

3,6,11,12,13,16,17,20 

Final  Number  of  GASP  Observables:    970 

Percentage  of  Data  Rejected: 

1.4% 

RMS  of  Residuals: 

4.947  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2693400.0858  -4346212.4518  3799864.7808 

Standard  Deviations  (X,Y,Z)  (meters): 

0.6707  0.5258  0.3626 

Geodetic  Coordinates  (monument): 

<D  =  36°48'09.6671     A  = 

=  238°12'46.8266    h  =  -29.842  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.5673  meters 
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Station  name: 

LOBOS3 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Precise  WGS  84  EF##91034 

Satellite  Clock  States: 

Precise  PC91034 

Date  of  Occupation: 

Day  37,  1991 

GPS  Receiver: 

TI  4100  BEPP/CORE  (versions  3.7/5.1) 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

10.8  Hours 

PRN  Numbers  Tracked: 

3,6,11,12,13,16,17,20 

Final  Number  of  GASP  Observables: 

1364 

Percentage  of  Data  Rejected: 

0.7% 

RMS  of  Residuals: 

4.956  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2693398.8698  -4346211.9262  3799864.2645 

Standard  Deviations  (X,Y,Z)  (meters) 

0.6641  0.4467  0.3459 

Geodetic  Coordinates  (monument): 

<D  =  36°48'09.6748    A  =  : 

238°12'46.8571     h  =  -31.022  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1.5673  meters 
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Station  name: 

LOBOS3 

Method: 

GPS  WGS  84  Absolute  Point  Positioning 

Software: 

STARPREP  (version  1.1),  GASP  (version  2.0) 

Ephemerides: 

Precise  WGS  84  EF##91034 

Satellite  Clock  States: 

Precise  PC91034 

Date  of  Occupation: 

Day  38,  1991 

GPS  Receiver: 

TI  4100  BEPP/CORE  (versions  3.7/5.1) 

Date  of  Point  Position  Estimation: 

June  1991 

Data  Collection  Span: 

10.8  Hours 

PRN  Numbers  Tracked: 

3,6,11,12,13,16,17,20 

Final  Number  of  GASP  Observables: 

1298 

Percentage  of  Data  Rejected: 

0.4% 

RMS  of  Residuals: 

5.160  cm 

WGS  84  Estimated  Station  Coordinates 

Monument  (X,Y,Z)  (meters): 

-2693398.3618  -4346212.7487  3799864.9918 

Standard  Deviations  (X,Y,Z)  (meters) 

0.6893  0.4706  0.3577 

Geodetic  Coordinates  (monument): 

d>  =  36°48'09.6858    A  =  : 

Z38°12'46.8920    h  -  -30.241  meters 

Antenna  height  (monument  to  electrical  center  of  antenna): 

1 .5673  meters 
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Attention:  Mr.  B.  Klepezyuski 

24.  Naval  Air  Warfare  Center 
Weapons  Division 

Pt.  Mogo,  CA  93042-5000 
Attention:  Mr.  B.  Cohenour, 
Mr.  R.  Smith 

25.  Naval  Oceanographic  Center 
Stennis  Space  Center,  MS  39922-5001 
Attention:  Mr.  P.  Taylor  (Code  HS) 

26.  Naval  Surface  Weapons  Center 
Dahlgren  Laboratory 
Dahlgren,  VA  22448 
Attention:  Mr.  R.  Hill, 

Mr.  A.  Evans,  Mr.  B.  Herman  (K10) 

27.  Navy  Space  Systems  Activity  (Code  80) 
P.O.  Box  92960 

Los  Angeles  Air  Force  Station 
Los  Angeles,  CA  90009-2960 
Attention:  CMDR.  J.  K.  McDermott 
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